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SUMMARY 


This  report  presents  the  initial  efforts  to  develop  an  improved 
method  for  the  prediction  of  performance  and  stress  characteristics 
of  a  single  VTOL-type  propeller -wing -nacelle  combination  operating 
in  various  flight  conditions. 

*  The  problem  of  determining  the  airloads  is  treated  by  re¬ 

placing  the  propeller  blade  by  a  bound  vortex  of  unknown  strength  and 
«  approximating  the  wake  by  a  mesh  of  trailing  and  shed  vortex  elements. 

Engineering  beam  theory  is  used  to  model  the  elastic  be¬ 
havior  of  a  twisted  rotating  blade. 

An  iterative  approach  is  used  to  solve  the  aeroelastic  equations 
on  a  digital  computer. 

Experimental  values  of  thrust,  torque,  and  blade  strains 
obtained  on  the  propellers  of  the  0.  6-scale  model  of  the  XC-142  are  pre¬ 
sented,  along  with  theoretical  results  for  corresponding  flight  configura¬ 
tions.  However,  because  of  the  insufficient  accuracy  of  the  experimental 
data  collected,  no  definite  evaluation  of  the  model  developed  is  made. 
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INTRODUCTION 


With  the  advent  of  VTOL  aircraft,  more  stringent  requirements 
were  placed  on  the  designer's  ability  to  predict  accurately  the  static 
thrust  and  torque  characteristics  of  propellers.  The  static  thrust  and 
torque  had  to  be  predicted  much  more  accurately  than  was  heretofore 
considered  acceptable.  Now,  the  accuracy  with  which  these  quantities 
are  predicted  may  well  determine  whether  the  aircraft  will  fulfill  its 
anticipated  role  as  a  VTOL  vehicle.  Since  the  entire  lifting  capability 
rests  on  the  propellers,  a  small  underestimate  of  the  thrust  might 
greatly  reduce  the  payload. 

In  addition  to  the  performance  prediction  problems,  oscillatory 
stresses  of  unexpectedly  large  amplitudes  were  observed  on  the 
propellers  of  such  VTOL  aircraft  during  wind-tunnel  and  full-scale 
tests.  In  particular, these  stresses  were  noted  in  the  transitional 
range  from  a  hover  to  an  axial  flight  configuration. 

Propeller  designers  suddenly  found  themselves  reexamining 
"tried  and  true"  propeller  performance  and  stress  theories  and  finding 
them  to  have,  at  best,  questionable  applicability.  No  proven  theoretical 
approach  was  available  for  the  skewed  flight  range  in  which  high 
oscillatory  stresses  were  observed.  Furthermore,  they  found  that 
much  of  the  experimental  data  upon  which  previous  correlation  with 
theory  had  been  based  were  not  sufficiently  accurate  to  aid  in  the 
development  of  more  accurate  theories  or  in  the  extension  of  existing 
ones. 


The  U.  S.  Army's  recognition  of  these  problems  and  Cornell 
Aeronautical  Laboratory's  (CAL)  previous  work  in  rotary-wing  air¬ 
craft  (as  exemplified  by  References  1,  2,  3,  5,  6,  8,  9,  10,  and  11) 
led  to  the  effort  reported  herein.  The  purpose  of  the  effort  was  to 
develop  an  improved  method  for  predicting  the  performance  and  stresses 
of  VTOL-type  propellers  through  their  entire  flight  range.  The  basic 
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method  developed  by  Cornell  Aeronautical  Laboratory  and  reported  in 
Reference  11  for  the  prediction  of  rotary-wing  airloads  in  hovering  and 
forward  flight  was  modified,  extended,  and  applied  to  the  problem  of 
VTOL-type  propellers. 

The  program  undertaken  was  twofold.  On  the  one  hand,  a 
theoretical  method  was  developed,  while  on  the  other,  experimental 
data  were  gathered  in  an  attempt  to  provide  a  consistent  set  of 
performance  and  stress  data  of  high  accuracy  for  correlation  with  theory. 

Throughout  this  report  the  term  performance  will  be 
synonymous  with  the  mean  forces  and  moments  acting  on  the  shaft; 
namely, thrust,  torque,  normal  and  side  force,  and  pitching  and  yaw¬ 
ing  moment. 
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THEORETICAL  PROGRAM 


STATEMENT  OF  THE  AEROELASTIC  PROBLEM 

Solution  of  the  aeroelastic  problem  requires  a  knowledge  of 
the  interaction  of  the  aerodynamic,  elastic,  and  inertial  forces  acting  on 
a  wing  section.  Once  these  interactions  are  known,  the  loads  and 
motions  of  the  section,  and  hence  the  total  wing  loads  and  motions,  may 
be  determined. The  problem  is  a  difficult  one  for  conventional  wings. 

For  rotary  wings,  i.  e.,  propeller  blades,  the  problem  is  compounded 
by  the  nature  of  the  wake  shed  from  the  blades,  which  now  takes  on  the 
geometry  of  a  distorted  helix;  by  the  rotation-induced  elastic  coupling 
between  various  blade  motions;  and  by  the  complicated  inertia  loads 
which  arise  when  motions  occur  in  a  rotating  frame. 

In  addition,  there  are  interference  effects  from  wings  and 
nacelles  which  lie  in  the  slipstream  and  which  may  be  quite  close  to  the 
propeller  plane.  To  define  the  interaction  of  these  forces,  the  following 
aerodynamic  and  structural  models  were  employed.  (The  discussion 
will  be  limited  to  steady-state  flight.  ) 

Aerodynamic  Model 

The  aerodynamic  model  accounts  for  the  effects  on  the 
blade  of  (1)  a  distorted  wake,  (2)  a  wing  and  nacelle  lodged  in  the 
propeller  slipstream,  (3)  forward  flight,  and  (4)  blade  deformations 
under  load. 


The  entire  model  rests  upon  periodicity  and  the  ability  .o 
predict  the  wake  position  with  sufficient  accuraf  y. 

Representation  of  the  Blade 

The  blade  is  represented  by  a  bound  vortex  of  strength, 
r  ,  placed  at  the  blade  quarterchord.  The  T  of  the  bound  vortex  varies 
both  with  radius  and  with  time.  Since  sueady-state  flight  conditions  only  ar 
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being  considered,  the  time  variation  is  periodic  of  period  T '  T'is 
the  time  for  one  rotor  revolution;  hence,  the  time  variation  of  r 
is  equivalent  to  the  azimuthal  variation.  Thus,  the  radial  and  azimuth¬ 
al  variations  of  T  are  given  by  cLr/dX  and  cLr/dLp{  X  is  the  radial 
position  and  V  is  the  azimuthal  angle),  respectively.  The  radial 
distribution  is  represented  by  breaking  the  bound  vortex  up  into  NR 
spanwise  segments  of  constant  r  .  The  azimuthal  variation  is 
similarly  accounted  for  by  allowing  each  radial  segment  to  assume  a 
new  value  of  r  as  it  is  stepped  through  HA  azimuthal  stations.  Thus, 
the  radial  and  azimuthal  variation  of  r  is  represented  by  (  AJR  *  a 7a) 
values  of  r  .  The  magnitude  of  r  is  determined  by  requiring  that  the 
chordwise  boundary  condition  of  zero  normal  velocity  at  the  three- 
quarterchord  of  the  blade  be  satisfied. 

Representation  of  the  Wake 

Because  of  the  variation  of  C  with  X  and  ^  ,  a 
continuous  sheet  of  vorticity  is  deposited  into  the  wake.  Because  of 
the  steady -state  flight  conditions,  the  variation  of  r  with  X  and  y  is 
repeated  for  each  and  every  blade  revolution.  The  wake  is  represented  in 
a  fashion  analogous  to  the  representation  of  the  bound  T  distribution. 

The  continuous  wake  is  modeled  by  a  mesh  of  "trailing"  and  "shed" 
vortex  elements.  The  "trailing"  vortex  element  is  deposited  in  the 
wake  from  the  ends  of  the  segmented  bound  vortex  and  is  equal  to  the 
difference  between  adjacent  bound  f's.  A  "shed"  vortex  is  deposited 
in  the  wake  from  each  bound  vortex  element  as  the  blade  moves  from 
one  azimuthal  position  to  the  next.  Its  strength  is  equal  to  the  differ¬ 
ence  between  bound  T 's  at  successive  azimuthal  positions.  Thus,  if 
a  given  trailing  vortex  is  followed  back  into  the  wake,  it  is  found  that 
its  strength  varies  in  a  stepwise  fashion.  It  is  constant  for  the  interval 
corresponding  to  the  increment  between  successive  azimuthal  positions 
of  the  blade.  Similarly,  consider  the  shed  vortex  deposited  in  the  wake 
at  some  azimuthal  position.  It  va  ies  in  stepwise  fashion  along  its 
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span.  The  variation  corresponds  to  the  spanwise  variation  of  the 
bound  vortex  which  existed  on  the  blade  when  the  blade  was  in  that 
azimuthal  position. 


Figure  I.  PICTORIAL  REPRESENTATION  OF  BLADE  WAKE  MESH 


Thus,  the  continuous  sheet  of  vorticity  deposited  in 
the  wake  has  been  replaced  by  a  mesh  of  si  d  and  trailing  vortex 
elements  whose  strengths  are  proportional  to  the  strength  of  the  bound 
vortex.  This  mesh  is  repeated  in  a  continuous  fashion  once  per  revolu¬ 
tion  of  the  blade  and  for  each  blade  in  the  propeller  set. 

Now,  consider  the  problem:  "How  is  the  wake 
positioned  relative  to  the  blades  which  generated  it?"  The  wake  position 
is  determined  by  the  interplay  of  such  effects  as  (1)  self-induced 
velocities  of  the  wake,  (2)  free-stream  velocity,  and  (3)  interference 
velocities  from  wings,  nacelles,  fuselages,  etc.  ,  in  and  in  the  proximity 
of  the  slipstream.  The  present  method  attempts  to  account  for  the 
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effects  of  (1)  wake  contraction  and  radial  and  axial  variation  of  the 
wake  velocity,  (2)  free -stream  velocity,  and  (3)  the  perturbing  effect  of 
a  wing  and  nacelle  completely  immersed  in  the  slipstream. 

The  wake  contraction  and  axial  and  radial  velocity 
variation  is  estimated  for  hovering  flight  by  using  a  system  of  finite-core, 
ring  vortices.  The  theory  was  developed  at  CAL  under  Office  of 
Naval  Research  sponsorship  and  is  reported  in  detail  in  Reference  1. 
Briefly,  the  approach  replaces  the  slipstream  by  an  axial  system  of 
discrete,  finite-core,  ring  vortices.  The  system  is  treated  as  a 
starting  problem  by  periodically  releasing  a  new  vortex  ring  in  the 
reference  plane  while  allowing  the  individual  rings  to  move  under 
mutual  influence  according  to  the  laws  of  vortex  dynamics.  After  a 
sufficient  number  of  rings  have  been  released,  a  stabilized  trajectory  is 
determined  down  which  all  successive  rings  move.  The  solution  was 
obtained  using  a  digital  approach,  and  the  program  output  yields  a  time 
history  of  the  coordinates  and  velocities  of  all  vortex  rings. 

This  model  was  applied  to  the  propeller  slipstream 
problem  in  which  the  wake  is  treated  as  a  mesh.  The  tip-trailing  vortex 
was  assumed  to  move  down  the  surface  generated  by  the  stabilized  ring 
vortex  trajectory.  This  was  done  by  requiring  that  the  tip  vortex  move 
with  the  velocities  with  which  the  rings  move.  Thus,  at  the  tip  of  the 
blade,  the  trailing  vortex  was  deposited  in  the  wake  with  the  velocities 
of  the  ring  at  the  instant  the  ring  was  deposited  in  the  reference  plane 
which  now  corresponds  to  the  prop  plane.  This  tip  vortex,  then,  was 
allowed  to  move  for  the  interval  of  time  corresponding  to  ^/VA  rotor 
revolutions.  At  the  end  of  this  period,  it  would  assume  new  values  of 
velocities  corresponding  to  the  velocities  which  the  ring  experienced 
at  that  position  in  the  wake,  and  so  on  down  the  wake.  Since  the  ring 
velocity  distribution  is  axially  symmetric,  tip  vortices  deposited  at 
various  azimuthal  stations  all  move  down  the  wake  with  the  same 
velocity  time  history  but  are  displaced  in  time  with  respect  to  one 
another  (see  Figure  2). 
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Figure  2.  PICTORIAL  REPRESENTATION  OF  TIP  VORTEX 
MOVING  DOWN  RING  VORTEX  TRAJECTORY 


The  result  is  a  tip  helix  of  va  ring  pitch  and  diameter. 

To  place  the  mesh  interior  to  the  tip  elements,  it  was 
assumed  that  (1)  all  interior  elements  move  with  the  same  axial  velocity 
as  the  tip  vortex  element,  and  (2)  the  radial  velocity  of  the  interior 
elements  is  proportional  to  the  product  of  the  radial  position  of  the 
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element  and  the  radial  tip  velocity. 


There  is  no  mechanism,  at  present,  to  allow  for  a 
tangential,  or  swirl,  component  of  velocity. 

This  wake  model  will  be  referred  to  as  the  RVM  (ring 

vortex  model). 


The  wake,  as  positioned  by  the  RVM,  is  used  as  a  base 
from  which  small  perturbations  due  to  the  free  stream,  wing,  and  nacelle 
are  allowed  in  the  case  of  low  forward  speeds.  These  perturbations 
are  included  as  incremental  velocities  at  each  control  point  in  the  wake 
(as  determined  from  the  RVM). 

Thus,  in  the  hovering  condition,  a  distorted  wake  is 
used  which  includes  the  effects  of  contraction,  radial  and  axial  velocity 
variations,  and  interference  effects  from  wing  and  nacelle.  For  the 
initial  stages  of  transition,  i.  e.  ,  low  forward  speeds,  the  free  stream 
is  included  as  an  additional  perturbation.  For  higher  forward  velocities, 
the  model  postulated  above  breaks  down, and  the  wake  is  treated  as  a 
skewed  helix  of  uniform  pitch,  with  none  of  the  effects  of  contraction  or 
of  the  wing  and  nacelle.  However,  if  a  wake  shape,  i.  e„  position,  is 
known  for  a  given  flight  condition,  it  can  be  included  in  the  program. 

For  details  of  the  mathematical  model,  see 
Appendix  I,  page  55. 

Representation  of  the  Wing  and  Nacelle 


The  interference  effects  of  the  wing  and  nacelle 
can  be  treated  in  terms  of,  first,  their  direct  effect  on  the  angle  of 
attack  at  the  propeller  and,  second,  their  effect  in  distorting  the  wake. 
The  direct  effect  is  always  included  when  applicable.  The  second 
effect  is  included  only  in  the  hovering  and  low  forward  speed  transition 
cases. 
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To  model  the  nacelle,  it  was  assumed  that  its  effects 
could  be  adequately  approximated  by  treating  it  as  a  Rankine  body 
(source-sink  combination).  The  Rankine  body  length  was  taken  equal 
to  the  nacelle  length,  and  the  Rankine  body  diameter  was  taken  such 
that  its  projected  frontal  area  equals  that  of  the  nacelle. 

The  wing  was  simply  treated  as  a  bound  vortex  placed 
at  the  wing  quarterchord.  For  details  of  the  mathematical  model  of 
the  wing  and  nacelle,  see  Appendix  II,  page  60. 

Development  of  the  Mathematical  Model  of  the 

Propeller  and  Wake 

The  equations  which  will  eventually  give  the  aero¬ 
dynamic  forces  required  for  the  aeroelastic  problem  will  now  be 
developed.  The  equations  will  be  formulated  in  terms  of  the  strength 
of  the  bound  vortex  r  .  Consider  an  element  of  a  blade  of  a  propeller 
in  skewed  flight  with  a  wing  and  nacelle  immersed  in  the  propeller 
slipstream  (see  Figure  3). 


PROPELLER,  WING,  AND  NACELLE 
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The  lift  per  unit  length  on  this  blade  element  is 


\ 


Jj  • /> VjO  (*> 

where 

/?  -  fluid  density 

Vj  -  local  total  velocity 

Pj  =  local  vortex  strength 

But  7;  can  also  be  written  as 

Jtj  --  bjfiV/cj.  (2) 

where 

c,  =  local  section  lift  coefficient 

J 

bj  =  local  blade  semichord 
Then  from  (1)  and  (2), 


0  *  >>j  vi  %  ■  (3) 

The  section  lift  coefficient,  c».  ,  is  a  nonlinear  function 

j 

of  the  effective  angle  of  attack,  ae-.  To  allow  for  this  nonlinearity,  as 
well  as  to  provide  a  mechanism  by  which  Reynolds  number  and  Mach 
number  effects  can  be  introduced,  the  c^.  curve  was  expanded  in  a  Taylor 
series  about  an  arbitrary  reference  angle  &0j-  ;  i.  e., 


Cj.  * 
*4 

dcA 

,  ,  /  d2Cj 

(<*<  -«*.)♦  -j  doe. 

It  is  now  required  that 

(<*e  ' 

J 


and,  hence,  terms  higher  than  (  0Ce  -  OcQ  ).  may  be  neglected  in 

J 

Equation  (4). 
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Thus,  Equation  (3)  becomes 


0  -  bj  Vj 


dc* 
C/  +  — — 
*•  dct 


(&e  ’  a*) 


Equation  (5)  is  regrouped  to  get 


where 


r*j  •  bJ  K 


Cj  -p- 

*  dct 


<*0 
(X,  J  J 


(5) 


(6) 


dCm 

a,  ® 

J  da 


a.. 


The  assumption  is  now  made  that 

sin  ao,  .  «  oc, .  . 

eJ 

This  approximation  leads  to  a  maximum  error  in  lift  of  approximately 
2%  when  o>.=  20°.  Since  as.. will  generally  be  much  less  than  20°, 

j  J 

it  was  felt  that  this  assumption  would  not  significantly  limit  the 
validity  of  the  theory. 

From  Figure  3,  it  is  recognized  that 

V  <X,  ~  V,  sin  Oc.  .  =  ir,  . 

From  lifting-line  theory,  Vj  must  be  the  velocity  perpendicular  to  the 
airfoil  chord  at  the  three-quarterchord  point.  Its  sense  must  be 
positive  up  to  give  increased  lift  for  increased  :  j  .  Equation  (6)  now 
becomes 


n  =  C 


*  bJ 


(?) 
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*  r^v  -  ' 


The  perpendicular  velocity  includes  contributions 
from  free  stream,  wake,  nacelle,  wing,  and  blade  elastic  motions;  thus, 

t*  =  [(X a  f  Vf  si n  T sin  y )  sin  6  -h  +  irw  +  v**  V’T]-  (8) 

where,  for  the  blade  element, 

X;  =  radial  distance  from  centerline  of  rotation  to  jth  blade 
element 

fi  -  rotor  speed 

dj  =  total  geometric  angle  including  built-in  twist,  jSj  ,  and 

elastic  motions,  <t>j ,  of  the  j**  blade  element 

hj  =  Vf  cosTcoaO  +  elastic  motions  perpendicular  to  the  airfoil 
chord  of  the  j** blade  element 

V*  -  velocity  due  to  wing  perpendicular  to  chord  of  j **  blade  element 

if*  -  velocity  due  to  nacelle  perpendicular  to  chord  of;<Ablade  element 

IT*  =  velocity  due  to  wake  perpendicular  to  chord  of  /*blade  element. 

Using  the  Biot-Savart  law  applied  to  the  wake 
mesh,  the  velocity  perpendicular  to  the  airfoil  chord  induced  by  the 
wake  can  be  written  as 


*  Z  °-ij 


(9) 


where 

-  the  geometric  portion  of  the  Biot-Savart  law  for  the  effect  of 
the  ith  element  at  the  jtH  point 

Y;  =  strength  of  the  itA  wake  element  of  vorticity. 


page  64. 


For  a  detailed  derivation  of  Equation  (9),  see  Appendix  III, 


Recall  that  the  / 's  of  the  wake  are  proportional  to  the 
r's  of  the  blade  bound  vortices  at  various  radial  and  azimuthal 
positions.  Hence,  when  Equation  (9)  is  substituted  into  Equation  (8) 
and  the  wake  terms  regrouped,  the  following  is  obtained  from  Equation  (7), 
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(10) 


QO 

[bJ  s  ^  -  *]  ri  +  b.ajYL  <Tijri  -  'tj 

i*  1 

where 

Jj  *  /^.  v-  bj  CLj  \^(XCL  +  Vf  sin  Z Sin  <1* )sin  6 -h  *■ 

and  CTjj  results  from  a  regrouping  of  the  <^j's. 

Equation  (10)  was  written  for  one  blade  element  at 
one  azimuthal  position.  Since  /j  changes  with  radius  and  azimuth, 
similar  equations  for  all  values  of  /j  can  be  written.  Since  /"}  is 
limited  to  (  NR  k  a /A  )  values,  (NR*NA  )  simultaneous  algebraic  equations 
in  /j  result. 


In  matrix  notation, 

[O  { r  f  -  W 


where 

is  a  square  matrix  of  order  (TTa^TIr)  whose  diagonal  elements 
have  the  form 

°Tj ' '  ] 


and  whose  off-diagonal  elements  have  the  form 


C7-.  • 


[Tj  is  a  column  matrix  of  the  unknown  /j 
[j]  is  a  column  matrix  of  the  known  quantities 


-  ro.  +  b-  <Li  | (X/1+  Vf  sin  Zstn  (p)  sin  9  -  h  ♦  ir*  +  ir  "J  J 
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At  this  point  in  the  development,  let  it  suffice  to 
say  that  everything  in  Equation  (11)  is  known,  except  the  /”j  's.  Hence, 
if  the  matrix  Equation  (11)  can  be  solved  for  the  Pj  's  ,  then  sufficient 
information  will  be  available  to  generate  the  required  aerodynamic 
forces  in  the  aeroelastic  problem. 

Summary  of  Assumptions 

To  develop  Equation  (11),  the  following  assumptions 

have  been  made: 

(1)  The  propeller -wing -nacelle  combination  is 
in  steady-state  flight. 

(2)  The  propeller  blade  can  be  represented 
aerodynamically  as  a  bound  vortex  placed  at 
the  quarterchord. 

(3)  The  radial  and  azimuthal  variation  of  r  can  be 
represented  by  a  finite  series  of  discrete  values. 

(4)  The  continuous  wake  can  be  represented  as 
a  mesh  of  concentrated  shed  and  trailing 
vortex  elements. 

(5)  sin  ae  ss  <xe  . 

(6)  The  wake  position  can  be  computed. 

(7)  The  wing  can  be  represented  by  a  concentrated 
vortex  at  its  quarterchord. 

(8)  The  nacelle  can  be  represented  as  a  Rankine  body. 

Structural  Model 

To  determine  the  inertia  and  elastic  forces  required 
for  the  aeroelastic  problem,  a  Lagrangian  approach  was  used. 

Following  the  development  of  Reference  7,  the  kinetic  and  potential 
energy  of  a  twisted  rotating  beam  was  written.  A  dissipation  term 
proportional  to  the  elastic  displacement  and  in  phase  with  the  elastic 
velocity  was  written  to  account  for  structural  damping  forces.  The 
elastic  displacements  were  then  expanded  in  terms  of  the  normal  modes 
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of  the  blades,  and  the  Lagrangian  was  applied.  Since  the  forcing  function 
of  the  system,  i.  e„  aerodynamic  load,  is  periodic  with  its  fundamental 
frequency  being /L  ,  the  steady- state  response  of  the  linearized  system 
must  be  periodic  of  fundamental  frequency  fl  .  Hence,  both  the  motions 
and  the  generalized  forces  were  expanded  in  a  Fourier  series  on  /l. 

Then,  by  the  uniqueness  theorem  for  trigonometric  series,  corresponding 
coefficients  must  be  equal.  The  result  is  a  set  of  (where  is  the 
number  of  modes  retained)  simultaneous  algebraic  equations  in  the 
Fourier  coefficients  (of  sine  and  cosine  components)  of  the  motions  for 
each  of  p  harmonics.  The  zeroth  harmonic  has  only  M  simultaneous 
equations,  since  there  are  no  corresponding  sine  terms  in  the  series. 

For  the  purpose  of  this  study,  six  degrees  of  freedom 
were  used— the  first  four  coupled  bending  modes  and  the  first  two 
uncoupled  torsion  modes.  A  detailed  development  of  these  equations 
and  the  results  are  presented  in  Appendix  IV,  page  69.  A  detailed 
developement  of  the  generalized  forces  is  presented  in  Appendix  V,  page  79. 

The  following  assumptions  were  made  in  defining  the 
structural  model: 

(1)  Simple  beam  theory  was  applicable. 

(2)  Airfoil  cross  sections  were  symmetrical  with  respect 
to  the  major  neutral  axis;  hence,  the  center  of  gravity 
lay  on  this  axis. 

(3)  Gravity  forces  were  negligible. 

(4)  Flight  conditions  were  steady  state. 

(5)  Damping  (structural)  was  proportional  to  the  elastic 
displacement  and  in  phase  with  the  elastic  velocity. 

(6)  Blade  elastic  deformations  we  e  small. 

(7)  The  blades  were  effectively  cantilevered,  and  the 
hub  had  infinite  impedance. 

All  the  forces  required  for  the  aeroelastic  problem  of 
a  propeller  in  transitional  steady-state  flight  have  now  been  defined. 

There  now  remains  the  problem  of  solving  the  resultant  set  of 
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simultaneous  equations  for  the  loads  and  motions. 


SOLUTION  OF  THE  AERO.  .ASTIC  PROBLEM 


Iterative  Solution 

An  iterative  approach  to  the  solution  was  adopted.  First, 
the  r  -  equations  are  solved  for  an  assumed  set  of  motions.  The 
resultant  loads  are  then  used  to  compute  the  generalized  forces  for 
the  equations  of  structural  deformation  of  the  blade.  The  equations 
of  blade  deformation  are  treated  as  a  forced  response  problem,  and 
new  motions  are  determined.  Since,  in  general,  the  new  motions 
do  not  agree  with  the  assumed  set,  the  f  -equations  must  be  resolved 
using  the  new  motions.  The  process  is  repeated  until  the  r  's  and 
blade  elastic  motions  from  the  ^iteration  agree  with  the  r  's  and  blade 
elastic  motions  from  the  (  n  +■  l)4*  iteration  to  within  some  specified 
tolerance.  At  this  point,  it  is  argued  that  a  consistent  set  of  loads 
and  motions  have  been  obtained,  since  no  change  in  either  occurs 
with  iteration.  Further,  since  the  problem  has  a  unique  solution, 
this  must  be  it. 

The  iterative  path  between  the  loads  (the  r  -equations) 
and  motions  (the  structural  equations)  is  interrupted  because  of  the 
nonlinear  character  of  the  curve.  Recall,  from  page  10,  Equation  (4), 
that  the  was  expanded  in  a  TayJ.or  series  about  some  known  reference 
angle  of  attack  cc0  .  So  long  as,  from  iteration  to  iteration,  oc0  <xe , 
to  within  some  amount  (*daf,  say  ),  the  same  &0  ,  CL,  and  are  used 

in  the  /’-equations  in  each  iterative  step.  If  <X0*  6  oc  <  oc  e  or  ac0-Aoc>Qce, 
a  new  oiQ  is  determined,  and  hence  a  new  a>  and  c \£  ,  for  use  in  the 

r  -equations  in  the  next  iterative  step.  This  is  accomplished  by 
dividing  the  c g  vs.  oc  curve  into  a  series  of  increments  in  CL  .  At 
the  center  of  each  is  evaluated  CC0  ,  a-,  and  c ^  .  Then,  as  each  new 

atg  is  determined,  it  is  checked  to  set  if  it  lies  in  the  range  corre¬ 
sponding  to  the  &0  used  to  compute  it  (i.  e. ,  does  <Xe  -ACCS  cc^  ^  <2^  ♦  A®). 

If  not,  it  is  placed  in  the  proper  range, and  the  iteration  is  continued.  The 
process  is  described  in  the  logic  diagram,  Figure  4. 
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Figure  *4.  LOGIC  DIAGRAM 


To  insure  convergence  of  the  iterative  scheme  employed 
to  solve  the  equations  of  motion, it  was  found  that  the  right-hand  side, 
i.  e.  the  generalized  force,  should  be  relatively  free  from  a  dependence 
upon  the  dependent  variables,  i.  e.,  the  elastic  motions.  In  the 
present  development  it  is  impossible  to  extract  this  dependence  of  the 
generalized  force  on  the  motions.  It  was  proposed  to  accomplish  this 
artificially  by  subtracting  from  both  sides  of  the  equations  of  blade 
deformation  the  contribution  from  an  approximation  to  the  quasi - 
steady  damping.  This  term  was  chosen  because  it  is  one  of  the 
largest  motion-dependent  contributors  to  the  generalized  force. 


The  term  to  be  subtracted  from  both  sides  of  the  equations 
of  blade  deformation  is  derived  in  Appendix  VII,  page  94.  The  treatment 
of  this  term  in  the  equations  of  blade  deformation  differs  depending  upon 
whether  we  consider  the  right-  or  the  left-hand  side  of  these  equations. 
On  the  right-hand  side,  the  motions  are  assumed  known  from  the 
previous  iterative  step.  On  the  left-hand  side,  the  deformations  are  the 
unknowns  to  be  used  in  the  succeeding  iterative  step.  For  example, 
consider  the  follov,  ng  representation  of  the  deformation  equation  at  the 
( n  +  1 )  iteration: 


mX 


n+1 


+  \cX 


n*  t 


+  CX 


n  +  i 


where  Fn  is  assumed  to  be  a  linear  function  of  the  elastic  motions 
determined  in  the  nih  iteration;  i.  e.  , 

Fn-  C°rC,*r,  "'Vo  ' 

Subtract  from  the  left-hand  side  and  CtX„  from  the  right-hand 

side  of  the  equation  to  obtain 

When  the  iterative  process  has  converged, 

-- 
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where 


6  £  -  predetermined  allowable  error. 

Then 


m, 


'^n+1  *  ^ ^ n-»-f  +  C^n+1  ~  > 

and  the  solution  is  effected. 


Performance  and  Stresses 

Once  a  converged  set  of  T's  and  motions  have  been 
obtained,  there  only  remains  the  proper  combination  of  these  quan¬ 
tities  to  obtain  the  forces  acting  on  the  propeller  shaft  and  the  blade 
stresses.  The  prop  thrust,  torque,  normal  force,  side  force, 
pitching  moment,  and  yawing  moment  are  then  determined  (see  Figure  5). 


Section  lift,  drag,  and  torsional  moment  (l-t  dj  ,  rr>4  respec¬ 
tively)  are  computed  at  each  radial  and  azimuthal  position  of  the  blade 
based  on  the  converged  loads  and  motions.  These  quantities  are  resolved 
into  components  with  respect  to  the  inertial  coordinate  system  (X0,  Y0  , 

ZJ.  At  each  azimuthal  position,  a  summation  of  the  spanwise  integration 
over  each  blade  yields  the  instantaneous  values  of  the  thrust,  torque,  etc. 
This  process  is  repeated  for  all  azimuthal  positions.  The  result  is  a  set 
of  NA  values  of  each  of  the  performance  quantities.  Hence,  with  A /A  equally 
spaced  azimuthal  positions,  the  performance  quantities  may  be  repre¬ 
sented  by  a  finite  Fourier  series  of  NA  terms,  the  zeroth  term  of  which 
yields  the  desired  mean  values  over  one  revolution: 


THRUST  = 


NA  r  a/C  [  NR 

T-^z[zzy^, 

k-t  Ljm' 


(12) 


NA  f  NC 

TotQue  = 


nr 


(13) 


NA  \TC 

NORMAL  FORCE  =  NF  -  +  = 

NA  k-.f  [«/*/ 


NR 


Y_  (f*j  c°s  v "  $ctj  +)&** 

J  --1 


4,fic 


(14) 
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SIDE  FOR.ce 


(15) 


HA 

a  5  F  -  ^  F 


hjc 

L 

4.S  / 


’a75 


y  (Fx.  sin  +  Fy-  coi 

j*f 


HA  MC 


PlTChilNG  MOMENT  -  PM^-=LY  <  V 

A/A f-  “ 


2^(^Xjeo»^  -  mj  sin  ^)^Xj 
As/  L^*'L^s/ 


(16) 


YAW  IMG  MOMENT  =  YM  -  =4  ^ 


NA 


RC 

1 

A»/ 


MR 


X"  ^jSinY^^cos^l IX, 


j *i 


4/J 


(17) 


where 


k  =  azimuthal  station  number 
=  blade  number 

/Vtf  =  total  number  of  azimuthal  stations 

tfC  -  number  of  blades 

Fx-  ~  '*&,(  dux/d  X)j  on  blade  a, 

fY;  -  ~  [dsin  ocj  +  JL cos  <4  on  blade  4. 

J4> 


[j cos  a i  -  JL  sin  &zj  -  on  blade  Ay 

2b*/>Vicm  ]+  E^lcos  &€  ■*■  dsin  &e]j  on  blade 
lj 

section  moment  coefficient  about  airfoil  quarterchord. 


A  detailed  development  of  the  stresses  at  any  point  in  the  blade  is 
presented  in  Appendix  VI,  page  86.  Consistent  with  the  structural  model 
adopted,  i.  e.,  simple  beam  theory,  only  the  longitudinal  and  shear 
stresses  are  developed  at  any  point  in  the  blade  cross  section. 

These  quantities  are  also  harmonically  analyzed. 
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The  Digital  Program 


The  entire  iterative  approach  to  the  solution  of  the 
aeroelastic  problem  has  been  programmed  in  FORTRAN  IV  for 
use  on  the  IBM  7044  digital  computer.  The  program  consists  of 
six  main  parts. 


Part  1 


This  program  positions  the  wake.  The  coordinates 
of  the  elemental  wake  vortex  end  points  are  computed.  These 
coordinates  are  computed, including  the  effects  of  ( 1 ) free-stream 
velocity,  (2)  wake  contraction,  (3)  wing  interference,  and  (4) 
nacelle  interference  in  any  combination  desired.  The  midpoints  of 
the  bound  vortex  elements  are  also  computed.  The  program  can 
handle  24  azimuthal  stations  and  10  radial  stations. 

Part  2 

From  the  coordinates  obtained  in  Part  1,  the 
effect  of  the  wake  vortices  of  unit  strength  on  the  blades  is  computed, 
i.  e.,  the  cr  of  each  individual  wake  element.  The  proper  combina¬ 
tions  are  made,  and  the  elements  which  contribute  to  the  [  o~  ]  matrix 
of  Equation  (11)  are  computed. 

Part  3 


Part  3  evaluates  all  the  mass  and  elastic  integrals 
necessary  for  computing  the  elements  of  the  structural  deformation 
equations.  It  also  writes  the  coefficients  of  the  least-square  poly¬ 
nomial  fit  to  the  two-dimensional  airfoil  section  data  used  to 
evaluate  ocg  ,  a,  ,  c*  ,  cmt  and  .  The  mass-elastic  integrals 
can  be  evaluated  at  a  maximum  of  30  spanwise  stations. 
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Part  4 


From  the  aerodynamic  data  of  Part  3  and  the  o~ 
elements  of  Part  2,  the  [  cr  ]  matrix  of  Equation  (11)  is  formed.  The 
set  of  simultaneous  equations  is  solved  by  using  the  Gauss-Siedel 
reduction  technique  (Reference  13).  The  technique  is  one  of 
iteration  until  the  set  of  r 's  from  one  iteration  agrees  with  the 
previous  set  to  within  some  prespecified  tolerance.  The  generalized 
force  is  then  computed  and  approximated  by  an  n order  polynomial 
in  X  •  This  is  required  by  the  fact  that  the  lift  is  computed  at  only  10 
spanwise  stations,  while  the  mass-elastic  properties  are  known  at  up 
to  30  stations. 


The  generalized  forces  are  harmonically  analyzed 
for  a  fundamental  period  of  one  propeller  revolution,  and  the  equa¬ 
tions  of  deformation  for  each  harmonic  are  solved  by  employing  the 
Gauss  elimination  process  (Reference  13). 

Overall  convergence  tests  on  ["  s  and  deformations 
are  made  by  comparing  values  obtained  in  the  iteration  to  those 
obtained  in  the  (n+  f)*h  iteration. 

If  the  r's  and  deformations  are  not  converged, 
and  ocI  are  recomputed  and  the  aerodynamic  tapes  are  searched  to 
determine  if  the  same  &0  is  to  be  used.  Control  is  returned  to  Part  3 
or  continues  to  Part  5,  depending  upon  the  results  of  the  convergence 
tests. 

Part  5 

Part  5  computes  a  time  history  of  the  propeller 
thrust,  torque,  side  and  normal  forces,  pitchin0  and  yawing  moments, 
and  stresses.  These  quantities  are  also  harmonically  analyzed.  The 
printed  output  from  the  program  includes  the  following  distributions 

given  as  functions  of  the  radial  and  azimuthal  position  in  the  disk;  (1) 
strength  of  bound  vorticity,  (2)  section  lift,  drag,  and  moment,  (3) 

section  displacement  (total),  (4)  section  effective  angle  of  attack, 
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(5)  section  induced  velocity,  and  (6)  section  stresses  —  longitudinal  and 
shear.  Also  printed  is  the  azimuthal  distribution  of  (7)  thrust,  (8) 
torque,  (9)  normal  and  side  forces,  and  (10)  pitching  and  yawirg  moments. 
The  harmonic  components  of  (1),  (2),  (3),  (6),  (7),  (8),  (9),  and  (10) 
are  also  printed. 
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EXPERIMENTAL  PROGRAM 


OBJECTIVE 


The  expressed  purpose  of  the  experimental  phase  of  this 
program  was  to  obtain  a  consistent  set  of  performance  (mean  thrust, 
torque)  and  stress  data  with  a  high  degree  of  accuracy.  These  data 
could  then  be  used  to  correlate  with  those  of  the  theoretical  phase 
of  this  program.  Of  the  available  existing  data  examined,  none  of 
the  performance  quantities  were  known  to  be  within  the  desired 
accuracy.  Furthermore,  no  set  of  stress  data  taken  in  conjunction 
with  performance  data  existed. 

At  the  time,  the  most  economical  means  for  obtaining  such  a 
set  of  data  seemed  to  lie  with  the  tests  of  the  0.  6-scale  model  of  the 
XC-142;  the  model  had  four  sets  of  three-bladed  propellers.  The 
tests  were  being  conducted  by  National  Aeronautics  and  Space 
Administration  -  Ames  in  the  40-ft  x  80-ft  wind  tunnel  for  the  purpose 
of  obtaining  individual  propeller  shaft  loads,  as  well  as  ship  total 
loads,  for  various  flight  configurations.  It  was  arranged  that  CAL, 
under  USAAVLABS  sponsorship,  would  be  allowed  to  instrument  one 
blade  of  each  of  two  prop  sets  in  order  that  bbde  strains  might  be 
measured  during  the  tests,  along  with  the  pr  seller  shaft  loads 
being  obtained  by  NASA.  Thus  a  consistent  set  of  strain  and  thrust 
and  torque  data  could  be  obtained,  recognizing  that  these  data  pertain 
only  for  a  specific  configuration  having  such  interference  effects 
as  wing,  wing  tip,  nacelle,  wing  flaps  and  slats,  propeller -propeller 
overlap,  and  fuselage.  Some  of  these  effects  could  have  been 
evaluated  by  performing  tests  with  selected  prop  sets  running.  For 
instance,  prop-prop  interference  could  be  eliminated  simply  by 
running  tests  with  only  one  prop  operative.  However,  due  to  test 
schedule  pressures,  no  such  configurations  were  tested. 
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MODEL,  APPARATUS,  AND  PROCEDURE 


NASA  was  responsible  for  all  aspects  of  the  tests  except 
gathering  of  the  strain  data,  which  was  CAL's  responsibility,  i,. 

Ames  was  to  collect  measured  values  of  prop  thrust,  torque,  normal 
force;  and  pitching  moment  by  using  a  rotating  balance  system  developed 
and  tested  by  NASA.  Since  both  the  performance  and  strain  data 
were  to  be  collected  in  a  rotating  frame  of  reference,  a  set  of  slip- 
rings  was  required.  Unfortunately,  the  number  of  rings  available 
was  insufficient  to  handle  both  the  strain  and  performance  data 
simultaneously,  and  a  switching  arrangement  to  be  mounted  on  the 
rotating  hub  had  to  be  devised.  So  long  as  the  flight  conditions,  i.  e., 
free-stream  velocity,  prop  rotor  speed,  etc.  ,  were  maintained 
constant  for  sufficient  time,  no  difficulty  was  anticipated  with  this 
sequenced  recording. 

Relays  to  provide  the  required  switching  capabilities  were 
tested  to  insure  their  operability  in  the  high  y  -field  to  which 
they  would  be  subjected.  These  tests  also  provided  the  assurance 
that  the  relays  would  not  interfere  with  the  electrical  signals  being 
transmitted  through  them  by  excessive  attenuation  or  by  contributing 
noticeable  amounts  of  noise  from  sources  such  as  relay  point  chatter. 

One  blade  on  each  of  two  prop  sets  was  instrumented  with 
two  rosette  strain  gages.  One  rosette  was  located  at  2.  67  feet  from 
the  centerline  of  rotation,  and  the  other  at  3.  67  feet  from  the  center- 
line.  Both  were  at  the  50%  chord.  These  locations  were  selected 
because  experimental  data  reported  in  Reference  2  indicated  that  the 
maximum  steady  bending  and  first  vibratory  bending  mode  strains 
would  be  maximum  at  50%  span  (  ~  2.  67  ft)  and  that  the  first 
vibratory  torsion  mode  strains  would  be  maximum  at  7  5%o  span 
(  3.  67  ft).  Since  the  blades  of  Reference  2  were  similar  to  the  blades 

on  the  0.  6 -scale  model  of  the  XC-142,  these  stations  were  selected. 
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Signal  wires  were  bonded  to  the  blade  trailing  edge  and  ter¬ 
minated  at  a  junction  plate  which  held  the  dummy  resistors  required  for 
the  gage  bridges  and  the  relays  required  for  switching.  From  the 

junction  plate,  signal  wires  were  attached  directly  to  the  platter-type 
slipring,  from  whence  they  entered  the  wind-tunnel  electrical  system. 
On  the  output  end  of  the  wind-tunnel  system,  parallel  recording  units 
were  provided  which  allowed  both  CAL  and  Ames  fo  record  the  rotating 
balance  signals.  Only  CAL  could  record  the  strain  signals. 

The  test  procedure  for  those  points  of  mutual  interest  to 
CAL  and  Ames  was  as  follows: 

(1)  Wind-tunnel  variables  (i.  e.  ,  tunnel  density,  velocity, 
total  force  measurements,  etc.  )  were  recorded  by 
Ames. 

(2)  Ames  recorded  the  signals  from  the  rotating  balance. 

(3)  The  signals  from  the  rotating  balance  were  then 
switched  to  and  recorded  by  the  CAL  units  (Arves 
units  being  shunted  out). 

(4)  The  relays  were  then  actuated.  This  allowed  one  set 
of  strain  data  to  be  recorded. 

(5)  The  relays  were  reset  to  a  new  position,  and  then 
the  second  set  of  strain  data  was  recorded. 

(6)  Control  was  then  returned  to  Ames. 

Just  prior  to  and  immediately  after  a  sequence  of  runs,  a 
resistance  calibration  check  was  made.  This  run  consisted  of 
switching  a  known  resistance  in  and  out  of  the  circuit  and  recording 
the  oscillograph  trace  deflection.  From  pretest  calibrations,  the 
load  per  inch  of  trace  deflection  could  then  be  d  .ermined. 

Also  recorded  at  the  beginning  and  end  of  each  run  were  the 
first  and  last  turnovers  of  the  propeller  with  the  wind  off. 

A  description  of  the  model  tested  is  given  in  detail  in 
Reference  4. 
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CONFIGURATIONS  TESTED 


The  prime  variables  of  interest  during  these  tests  were  wing 
angle  (  oew  )  with  respect  to  the  fuselage  centerline,  fuselage  angle 
(  otf  )  with  respect  to  the  free  stream,  prop  speed  (j"L),  tunnel  velocity 
(  *^),  and  prop  thrust  and  torque.  Also  of  interest  were  the  effects 
of  yaw  angle  ( <xr  ),  flap  angle  (ocp  ),  and  leading-edge  slats. 

A  configuration  run  was  determined  by  the  wing  angle.  For 
a  given  run,  subcases  were  determined  by  various  combinations 
of  the  above  variables.  A  configuration  was  referred  to  as  "clean" 
if  there  were  no  leading-edge  slats  and  <Xp  =  0.  The  configuration 
was  "dirty"  if  both  leading-edge  slats  and  flap  angle  were  employed. 

The  configurations  for  which  both  strain  and  performance 
data  were  collected  are: 


Clean  Configurations  (  5  0,  ocP  s  0)  (i.  e. ,  without 

flaps  and  slats) 


a*  =  40° 


c xF  *  10° 

OCp  »  0° 

OCp  =  -to*. 


Vf  =  39.8  FT/sec  i  n*  126*  r.p.m. 


oc *  -  60 


ocF  r  0°  i  Vf  -  33. 7  rT/5ec  J  n  «  1330  r.  p.  m. 


0CW=  90° 

OCF  =  o';  Vf  =  O  ;  fl  »  133  0  r.p.m. 
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and  slats) 


Dirty  Configurations  (dfr*0,  &p  =  60  )  (i.  e.  ,  with  flaps 


oc,w  = 


A 0 ° 

OCp 

<XF 

OCp 


10 •  1 


10°  ) 


Vf  r  39.  8  ft/sec;  »  1264  r. 
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DISCUSSION  OF  RESULTS 


THEORETICAL  RESULTS 

The  digital  program  discussed  on  page  22  was  used  to  obtain 
theoretical  results  for  seven  flight  cases  of  the  0.  6-scale  model  of  the 
XC-142.  The  mass-elastic  data  for  the  model  (i.  e.  ,  the  distributions 
of  mass,  stiffness,  inertia,  area,  etc.  )  were  obtained  from  the 
Hamilton  Standard  Company,  which  built  the  model  blades.  Also 
received  from  Hamilton  Standard  Company  and  incorporated  into  the 
program  were  two-dimensional  airfoil  section  characteristics  for 
the  airfoil  sections  used  in  the  model  blades. 

Five  of  the  cases  analyzed  correspond  to  the  clean  configura¬ 
tions  for  which  experimental  data  were  collected  (see  section  titled 
Configurations  Tested,  page  28-)  The  T  =  40°  case  was  computed  a  second 
time  to  show  the  effect  of  changing  the  wing  load  with  all  other 
parameters  held  fixed.  The  X  -  90°  (hovering)  was  computed  twice 
to  indicate  the  effect  of  ground  proximity  on  the  0.  6-scale  model  of 
the  XC-142.  The  time  histories  of  thrust,  torque,  normal  and  side 
forces,  pitching,  and  yawing  moments  are  presented  for  ^ach  case 
in  Figures  6  through  10. 

Note  that,  in  Figures  6  through  10,  the  mean  values  of  the 
thrust,  torque,  normal  force,  etc.  ,  are  printed  at  the  top  of  each 
graph  and  the  oscillatory  components  are  plotted  about  zero  mean. 

The  dominant  frequency  contents  of  the  oscillatory  values  of 
the  thrust  and  torque  for  all  cases  except  T  =  90°  are  3P  (i.  e.,  a 
frequency  equal  to  3  times  the  prop  speed),  6P,  9P,  etc.  The  peak- 
to-peak  magnitude  of  the  3P  component  of  the  thrust  is  approximately 
5%  of  the  mean  thrust;  the  magnitude  of  the  6P  component  is 
approximately  10%  of  the  3P  component.  The  magnitudes  of  the  IP, 

2P,  4P,  5P,  7P,  etc.  ,  harmonics  are  of  the  order  of  1  x  10  ^  lb. 

Since  the  only  harmonics  of  the  performance  quantities  which  can 
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be  transmitted  to  the  shaft  during  steady- state  flight  must  be  integer 
multiples  of  the  number  of  blades,  the  magnitudes  of  these  latter 
harmonic  components  become  a  measure  of  the  accuracy  of  the 
numerical  method  employed.  Torque,  normal  and  side  forces,  ;  ’ 

pitching  and  yawing  moments  exhibit  the  same  characteristics  as  the 
thrust  of  dominant  3P,  6P,  and  9P  harmonic s,  with  all  other  harmonics 
negligibly  small;  however,  for  the  normal  and  side  forces,  and  pitching 
and  yawing  moments,  the  phasings  of  the  significant  harmonics  are 
such  as  to  cause  their  time  (  if/  )  histories  to  appear  less  sinusoidal. 
Note  also  the  relatively  large  mean  and  oscillatory  pitching  and  yawing 
moments  predicted  in  all  cases. 


A  comparison  of  Figures  7  and  8  indicates  the  effect  of 
changing  the  wing  load.  With  all  other  parameters  held  fixed,  the 

2 

strength  of  the  wing  bound  vortex  {  rw  )  was  decreased  from  400  ft  /  sec 

(Figure  7;  this  case  corresponds  to  the  exp^  ^imental  configuration) 

2 

to  298  ft  /sec  (Figure  8),  a  reduction  of  2  /o.  Almost  no  change  in 
the  mean  thrust  or  torque  was  noted.  The  peak-to-peak  amplitudes 
of  the  time-dependent  portions  of  the  thrust  and  torque  were  reduced 
by  approximately  25%.  The  mean  values  of  the  remaining  performance 
quantities  (side  force,  pitching  moment,  etc.  )  were  reduced  by 
approximately  11%,  while  the  oscillatory  portions  of  these  quantities 
were  reduced  by  approximately  21%.  The  apparent  one-to-one 
correspondence  between  and  the  peak-to-peak  amplitude  of  the 
performance  quantities  is  believed  to  be  coincidental. 

Two  hovering  (  r  -  90°,  -  0)  cases  were  computed.  One 

case  was  intended  to  simulate  the  model  sitting  on  the  ground; 
hence,  the  prop  plane  was  approximately  one  ro.or  diameter  above 
the  ground;  the  other  case  was  computed  for  the  same  parameters 
except  with  the  prop  plane  well  out  of  ground  effect.  It  was  found  that 
the  effect  of  the  ground  was  slightly  detrimental.  The  table  below 
presents  the  thrust  and  torque  values  for  the  two  cases. 
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Thrust  (lb) 

Torque  (ft -lb) 

Hovering  In-Ground  Effect 

1176 

639 

Hovering  Out-  of-  Ground  Effect 

121  3 

642 

While  the  torque  remained  essentially  constant,  the  thrust 
decreased  by  approximately  3%, due  to  the  effect  of  ground.  Examina¬ 
tion  of  the  RVM  in  the  two  cases  revealed  that  for  the  near 
wake  (i.  e.  t  within,  say,  a  quarter  prop  radius)  the  shape  and  spacing 
of  the  wake  were  essentially  the  same.  Farther  down  the  wake,  how¬ 
ever,  the  shape  and  spacing  were  no  longer  the  same.  At  approxi¬ 
mately  one  prop  radius,  the  wake  boundary  for  the  in-ground  effect 
case  increased  in  diameter,  while  that  for  the  out-of-ground  effect 
case  continued  to  contract.  It  is  believed  that  for  the  case  in  hand 
the  small  increase  in  thrust  and  torque  obtained  as  the  height  above 
ground  increased  was  due  to  the  fact  that  the  far  wake  had  contracted 

more.  That  the  increase  was  small  is  attributed  to  the  fact  that  the 

❖ 

near  wakes  are  almost  identical  in  the  two  cases. 


Note  that  the  effect  of  the  ground  can  be  either  beneficial  or 
detrimental,  depending  upon  (l)the  degree  of  wake  distortion,  (2) 
the  wake  spacing,  and  (3)  the  load  distribution  on  the  blade.  For 
example,  if  the  blade  load  distribution  is  such  that  a  major  portion 
of  the  load  is  carried  on  the  most  outboard  blade  sections  (say 
from  90%  span  out),  then  contraction  of  the  wake  will  result  in  a 
decrease  in  load.  That  is  provided  that  the  wake  contraction 
does  not  move  inboard  of  the  90%  span.  If  the  contracted  wake 
moves  inboard  of  the  90%  span,  the  wake  effect  is  such  as  to 
increase  the  load. 

An  additional  contributing  factor  which  was  not  included  in  the  above 
analysis  is  the  influence  of  the  ground  image  system  on  the  induced 
velocity  at  the  rotor  plane. 
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Presented  in  Figure  11  are  the  spanwise  variations  of  the 
harmonics  (zeroth  to  the  fourth  only)  of  the  theoretically  computed 
longitudinal  stress  for  T  -  50°.  The  plot  is  representative  of  the 
stress  distributions  computed  for  all  the  cases. 

EXPERIMENTAL  RESULTS 

Considerable  difficulty  was  encountered  in  trying  to  reduce 
the  experimental  data.  Even  before  the  wind-tunnel  program  was 
completed,  peculiarities  were  noted  in  the  data  from  nacelle  number  1 
(the  outboard  nacelle  on  the  pilot's  left).  It  was  subsequently  found 
th  t  an  intermittent  short  had  developed  in  one  of  the  gages  of  the 
balance  &  stem.  Since  all  the  gages  on  nacelle  number  1  were 
powered  from  a  common  supply,  it  was  felt  that  the  signals  of  all 
gages,  both  from  the  balance  system  and  the  strain  gages,  would  be 
affected  by  the  intermittent  short.  Hence,  none  of  the  data  collected 
on  nacelle  1  wer>  .  uuced. 

The  raw  thrust  and  torque  and  strain  data  collected  from 
nacelle  number  2  (inboard  nacelle  on  pilot's  left)  were  examined.  It 
was  noted  that  these  were  not  periodic  over  the  time  interval  during 
which  records  were  taken.  This  nonperiodicity  was  much  more 
pronounced  on  the  strain  records  than  on  the  thrust  and  torque  records. 

It  was  determined  that  the  nonperiodic  loadings  were  due  to  the  fact 
that  the  two  adjacent  propellers  were  not  synchronized;  i.  e.,  they  did 
not  have  the  same  prop  speed.  Thus,  the  propellers  experienced  a 
slightly  different  wake  geometry  for  every  revolution.  Since  the  difference 
in  speed  was  small,  the  effect  of  the  wake  changes  was  not  as  notice¬ 
able  in  the  thrust  and  torque  data  as  in  the  strain  data,  which  showed 
variations  from  one  revolution  to  the  next  in  some  cases. 

Harmonic  analysis  of  the  thrust  and  torque  data  revealed 
relatively  large  IP  and  2P  frequency  components  (on  the  order  of  10% 
of  the  mean  values  for  some  cases).  These  components  could  not  have 
existed  if  the  thrust  and  torque  had  been  periodic. 
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Because  of  nonperiodicity  of  the  loads,  the  sequencing  procedure 
(see  section  titled  Model,  Apparatus,  and  Procedure,  page  26)  used 
in  collecting  the  test  data,  in  particular  the  strains,  becomes  highly 
questionable.  Reduction  of  some  of  the  strain  data  was  performed 
based  upon  certain  restrictive  assumptions  as  discussed  in 
Appendix  VIII,  page  97.  Since  the  thrust  and  torque  exhibited  the  least 
effect  of  the  nonperiodicity,  it  was  reasoned  that,  at  least,  the  means 
of  these  quantities  could  be  obtained  from  the  data. 

Reduction  of  the  thrust  and  torque  data  posed  another  problem 
in  addition  to  the  ones  concerned  with  nonperiodicity.  Certain 
anomalies  were  noted  in  the  rotating  balance  data  collected  during  the 
test  which  gave  rise  to  questions  concerning  the  accuracy  of  the 
calibration  of  the  balance  unit  and,ii  t  ticular,  to  the  amount  of 
crosstalk  between  the  various  components. 

The  omalies  were  observed  in  the  th.ust  and  torque  data 
taken  during  the  very  slow  rotation  of  the  propellers  prior  to  bringing 
the  wind  tunnel  up  to  speed.  Because  of  the  low  prop  r.  p.  m.  ,  it  was 
expected  that  the  inertia  and  aerodynamic  loads  would  be  negligibly 
small.  Since  the  weight  components  had  been  electrically  nulled  at 
zero  r.  p.  m.  ,  no  changes  from  the  zero  r.  p.  m.  positions  in  the  thrust 
and  torque  signals  were  expected  at  these  low  r.  p.  m.'s.  However,  a 
mean  shift,  plus  a  IP  oscillatory  signal  in  thrust  and  torque,  was 
observed.  The  origin  of  this  signal  is  unknown,  as  is  its  behavior  with 
the  variables  of  a  ,  Z  .  etc. 

Because  of  the  errors  associated  with  (1)  data  recording  and 
reduction,  (2)  the  nonperiodicity  of  the  loads,  and  (3)  the  above  noted 
anomalies  in  the  experimental  thrust  and  torque,  the  data  presented  in 
this  report  for  the  thrust  and  torque  are  believed  to  be  good  only  to 
within  ^  10%. 

Presented  in  Figures  12  through  17  are  the  mean  and  oscilla¬ 
tory  values  of  the  harmonically  analyzed  experimental  thrust  and  torque. 
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The  dominant  nonzero  harmonic  of  the  thrust  for  all  cases  including 
7  -  90°  is  the  3rd.  It  has  a  peak-to-peak  value  of  approximately  10% 
of  the  mean.  Note  the  large  IP  frequency  component  in  most  of  the 
configurations.  The  oscillatory  portion  of  the  torque  contains  large 
components  at  frequencies  from  IP  through  6P  in  almost  all  cases. 

A  comparison  of  Figures  14  and  15  shows  the  effects  of  adding 
slats  and  60  of  wing  flaps.  Figure  14  is  the  clean  configuration,  i.  e., 
without  slats  and  flaps,  while  Figure  15  is  the  dirty  configuration,  with 
slats  and  flaps.  The  mean  and  oscillatory  components  of  the  thrust 
show  almost  no  effect  due  to  slats  and  flaps.  However,  the  mean  torque 
was  reduce  y  approximately  10%  with  the  addition  of  flaps  and  slats, 
while  the  dominant  frequency  component  shifted  from  3P  to  6P. 

Figure  17  presents  the  thrust  and  torque  for  the  hovering  case. 
Note  the  relatively  large  oscillatory  components  of  thrust  and  torque. 
The  peak-to-peak  oscillatory  thrust  is  approximately  10%  of  the  mean 
thrust.  The  peak-to-peak  oscillatory  torque  is  approximately  7%  of 
the  mean  torque.  These  oscillatory  loads  are  believed  to  be  due  to 
the  interference  effects  of  prop-prop  overlap,  the  wing,  nacelle,  and 
fuselage. 


Because  of  the  nature  of  the  rotating  balance  and  the  fact 
that  the  prop  loads  were  not  periodic,  the  total  normal  force  could 
be  measured  only  at  y  =  0  and  180°,  while  the  total  side  force  could 
be  measured  only  at  <f>  =  90°  and  270°.  Hence,  no  information 
concerning  the  nature  of  the  mean  or  oscillatory  normal  forces  or 
side  forces  could  be  obtained  from  the  measured  data.  The  same 
comments  apply  to  the  pitching  and  yawing  moments. 

COMPARISON  OF  THE  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

Comparison  of  the  mean  theoretical  and  experimental  thrust 
and  torque  are  oresented  in  Figure  18.  The  experimental  data  values 
are  shown  with  a  band  indicating  the  estimated  error  of  *  10%. 
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In  the  range  30°  £  t  —  50°  the  theoretical  thrust  values  lie 
within  the  estimated  accuracy  of  the  experimental  thrust  data.  At 
7"  =  60°  and  T  =  90°  the  theoretical  thrust  is  slightly  greater  than 
the  experimental  thrust. 

The  theoretically  predicted  torque  was  greater  than  the 
experimentally  measured  torque  for  all  configurations  tested. 

No  comparison  between  the  oscillatory  loads  could  be 
made  because  of  the  nonperiodicity  of  the  experimental  thrust  and 
torque. 


Since  the  experimental  data  are  not  sufficiently  accurate 
and  include  interference  effects  (prop-prop  overlap  and  fuselage) 
not  included  in  the  present  mathematical  model,  it  was  not  possible 
to  make  sufficient  comparisons  for  a  definitive  evaluation  of  the 
theoretical  method  developed. 


36 


Figure  6.  MEAN  AND  INSTANTANEOUS  VALUES  OF  THEORETICAL  PERFORMANCE 
FOR  V*  =  39.8  ft  /sec  ;  Z  =  30*;  SI  =  1264  r.p.m.; 

rw  =  310  ft  2/sec 
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Figure  7.  MEAN  AND  INSTANTANEOUS  VALUES  OF  THEORETICAL  PERFORMANCE 
FOR  /f  =  39.8  ft  /sec  ;  z  =  40*;  si  =  1264  r.p.m.; 

-  400  ft  2/sec 
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Figure  8.  MEAN  AND  INSTANTANEOUS  VALUES  OF  THEORETICAL  PERFORMANCE 
FOR  Vf  -  39.8  ft  / sec  ;  r  =  40*;  fl  =  1 26 >  r.p.m.; 

/>  =  298  ft  2/sec 
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Figure  9.  MEAN  AND  INSTANTANEOUS  VALUES  OF  THEORETICAL  PERFORMANCE 
FOR  Vp  =  39.8  ft  /sec  ;  r  =  50°;  n  -  1 264  r.p.m.; 

r„  =  479  ft  2/sec 
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Figure  10.  MEAN  AND  INSTANTANEOUS  VALUES  OF  THEORETICAL  PERFORMANCE 
FOR  Vr  -  33.7  ft  /sec  ;  r  =  60*;  ft  -  1380  r.p.m.; 

rw  =  457  ft  2 /sec 
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Figure  II.  THEORETICAL  HARMONICS  OF  LONGITUDINAL  STRESS 
VS.  RADIAL  POSITION  FOR  T  =  50° 
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Figure  12.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPcRIMENTAL  THRUST  AND 

TORQUE  DATA  FOR  Vf  =  39.8  ft  / sec  ;  ?  =  30°;  -ft  =  1264  r.p.m 
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Figure  13.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  THRUST  AND 

TORQUE  DATA  FOR  Vf  =  39.8  ft  /sec  ;  r  =  W;  A  =  1264  r.p.m. 
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Figure  14.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  THRUST  AND 

TORQUE  DATA  FOR  Vp  =  39.8  ft  / sec  ;  r  =  50°;  Cl  =  1264  r.p.m 
WITHOUT  FLAPS  AND  SLATS 


THRUST,  1b  TORQUE,  ft  -lb 


Figure  15.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  THRUST  AND 

TORQUE  DATA  FOR  Vf  =  39.8  ft  /sec  ;  Z  -  50°;  il  =  1264  r.p.m.; 
WITH  FLAPS  AND  SLATS 


46 


THRUST,  1b  TORQUE,  ft  -lb 


0  80  160  2 >10  320 

V' ,  AZIMUTH  ANGLE,  deg. 

Figure  16.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  THRUST  AND 

TORQUE  DATA  FOR  =  33.7  ft  /sec  ;  Z  =  60';  a  =  1380  r.p.m. 
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Figure  17.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  THRUST  AND 

TORQUE  DATA  FOR  Vp  =  0  ft  /sec  ;  f  -  90*;  J1  =  1380  r.p.m. 
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Figure  18.  COMPARISON  OF  MEAN  THEORETICAL  AND  EXPERIMENTAL  THRUST  AND 
TORQUE  VS.  ANGLE  OF  ATTACK  OF  THRUST  AXIS 
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CONCLUSIONS 


Based  on  the  results  obtained  in  this  study,  the  following 
conclusions  were  drawn. 

It  has  been  shown  during  the  course  of  this  study  that  the 
approach  employed  to  define  and  formulate  the  aeroelastic  problem 
associated  with  VTOL-type  propellers  operating  in  all  flight  regimes 
from  hovering  (in-  and  out-of-ground  effect)  through  transition 
and  into  axial  flight  is  a  feasible  one.  Furthermore,  the  digital 
program  developed  to  implement  solution  of  the  resultant  aeroelastic 
equations  converges  to  a  unique  solution  in  a  reasonable  length  of  time 
(approximately  0.  75  hour  on  the  IBM  7044). 

The  experimental  thrust  and  torque  data  collected  during  the 
NASA-Ames  wind-tunnel  test  program  are  believed  to  be  good  only  to 
±  10%, which  is  not  considered  sufficiently  accurate  to  provide  the 
desired  evaluation  of  the  model  developed  herein.  Although  it  was 
possible  to  collect  the  thrust  and  torque  data,  certain  other  needed 
information,  such  as  the  mean  and  oscillatory  normal  and  side  forces  as 
well  as  pitching  and  yawing  moments,  was  not  obtained  because  the 
propeller  load  measurement  system  could  not  properly  account  for 
the  nonperiodic  character  of  the  data.  As  a  result,  the  rotating 
balance  unit  employed  in  these  tests  yielded  data  which  could  not  be 
converted  into  the  desired  inplane  forces  and  moments.  Furthermore, 
it  was  concluded  that  the  measured  blade  strains, which  were  also 
nonperiodic, could  not  be  converted  to  stresses  without  certain 
assumptions.  These  assumptions  were  highly  questionable  in  light 
of  the  nonperiodicity  of  the  strain  measurements;  hence,  the 
stresses  obtained  were  questionable. 
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RECOMMENDATIONS 


CAL  is  currently  working  on  a  modification  to  the  contract 
under  which  this  report  was  written.  The  purpose  of  this  modification 
is  to  obtain  accurate  experimental  thrust,  torque,  and  blade  strains 
and  stresses  on  an  isolated  propeller  ( i,  e.  ,  no  interference  effects 
from  other  props,  wings,  etc.  ).  The  data  will  be  used  to  correlate 
with  those  of  the  theoretical  approach  developed  herein.  This  effort 
is  expected  to  satisfy  the  foremost  recommendation,  namely,  that  of 
obtaining  accurate  experimental  data  for  an  isolated  propeller. 

It  is  further  recommended  that  (1)  the  computational  program 
be  employed  to  provide  quantitative  answers  to  the  effects  of  wing, 
wake,  forward  velocity,  and  nacelle,  (2)  the  present  program  be 
expanded  to  include  the  effects  of  prop-prop  overlap  and  fuselage,  and 
to  investigate  whether  the  lifting -line  model  of  the  blade  and  wing 
should  be  abandoned  in  favor  of  a  model  which  accounts  for  a  chord- 
wise  distribution  of  vorticity,  and  (3)  theoretical  and  experimental 
programs  be  undertaken  with  the  objective  of  improving  the  wake  model. 
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APPENDIX  I 


POSITIONING  THE  WAKE 

HIGH  AXIAL  AND  TRANSITIONAL  FORWARD  VELOCITIES 

For  the  case  where  the  wake  is  carried  back  from  the  prop 
plane  at  high  velocity  in  both  the  axial  and  in-plane  direction,  the  wake 
is  assumed  to  be  a  helix  of  constant  pitch  having  no  contraction.  The 
interference  effects  of  the  wing  and  nacelle  on  the  wake  shape  are 
neglected.  The  coordinates  of  any  point  in  the  wake  can  then  be 
determined  in  a  straightforward  fashion.  Consider  a  blade  in  some 
azimuthal  position  trailing  a  mesh  of  vorticity.  Let  KA  denote  the 
azimuthal  position  of  the  blade  (KA  -  i,  2,  ...  ffiA  ),  *  denote  radial 
position  of  a  blade  element  (  fc  =  1 ,  2,  ...  NR),  L  denote  the  azimuthal 
position  in  the  wake  (i  =  1  (at  the  blade),  2,  .  .  .  WAtW  where  IVis  the 
number  of  revolution  of  wake),  and  j  denote  the  radial  position  in  the 
wake  ( j  =1,2,  ...  /VR). 


Note:  At  the  end  of  each  appendix  is  a  list  of  symbols  which  applies  to 
that  appendix. 


Then  the  coordinates  of  any  point  in  the  wake,  for  a  given  rotor 
position,  are  given  by 


rj  cos  (jKA-i)A^)  +  V{  sin  T ^ “ 0 


(18) 


%  '  rj  3<»({KA-i}M)  (19) 

t;j  =■  \W  -  \/f  cos  0  (20> 


where 

fj  =  radial  distance  from  centerline  of  rotation  to  j*h  point  in 
wake 

W  =  inflow  velocity  (constant  over  disk)  . 

HOVERING  AND  LOW  TRANSITIONAL  VELOCITIES 

For  the  cases  of  hovering  and  low  transitional  velocities,  the 
effects  of  wake  contraction,  axial  and  radial  velocity  variation,  wing, 
nacelle,  and  forward  velocity  are  accounted  for.  The  effects  of  wing, 
nacelle,  and  forward  velocity  are  treated  as  perturbations  from  the 
wake  geometry  obtained  from  the  ring  vortex  model  (RVM).  The 
ring  vortex  model  represents  the  rotor  wake  by  a  series  of  discrete 
ring  vortices  which  are  deposited  in  the  wake  one  at  a  time  at  specified 
intervals  and  then  allowed  to  move  under  their  own  mutual  interactions. 
After  a  sufficient  number  of  rings  have  been  deposited  in  the  wake,  a 
stabilized  trajectory  is  obtained  down  which  all  subsequent  rings 
pass.  This  trajectory  is  used  to  represent  the  wake  trajectory  of  an 
equivalent  rotor.  The  RVM  yields  a  time  history  of  the  radial  and 
axial  components  of  displacement  and  velocity  of  each  ring  released 
into  the  wake.  This  information  is  used  to  position  the  propeller 
wake. 
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In  accordance  with  Figure  2  and  Figure  19.  the  coordinates  of  the  tip 
vortex  are  given  by 

Su,  *  if  [vfr(jf-)‘°s({KA +Rcos  (21) 

5 =  1 

d-t)  r  .  * 

(22) 

5«/  L  J 

«. [■»?(¥;] 

S  =  t 

where 

* 

n 

w  =  radial  velocity  of  the  ring  at  time  6 

l Vg*  =  axial  velocity  of  the  ring  at  time  6 

R  =  distance  from  centerline  of  rotation  to  blade  tip. 

To  position  the  wake  mesh  interior  to  the  tips,  it  is  assumed  that  the 
inboard  elements  of  vorticity  are  translated  with  the  same  axial 
velocity  as  the  tip,  while  the  radial  velocity  is  taken  as  proportional 
to  (-—  )  Vqj  .  Thus  Equations  (21)  through  (23)  become,  for  any  point 
in  the  wake, 


?ij‘Z  \$)v/r(^)c°s({M-i!M)j  +  rjcos({KA-i}M)  (24) 

<5*  / 

(*-/) 

nH  •L  [(^)^r(^-)scn([KA-i}^) ]  + rj  sln({KA-t]At)  (25) 

6 « 1  L 

■'?  [».?  &)] .  »> 

S  •  t 


Now  if  we  evaluate  the  perturbation  velocities  due  to  the  wing  and 

nacelle  at  the  points  (  ,  rj-j ,  £-j  )  and  denote  them  as 

Note:  Recall  that  for  the  steady-state  conditions,  time  and  azimuthal 
position  are  synonymous. 
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radial  velocity  at  j**  wake  point  at  time  S  due  to  the  nacelle 

axial  velocity  at  j*h  wake  point  at  time  8  due  to  the  nacelle 

X  component  of  velocity  at  j**1  wake  point  at  time  6  due  to  the  wing 

Y  component  of  velocity  at  j^wake  point  at  time  8  due  to  the  wing 

Z  component  of  velocity  at  j**  wake  point  at  time  8  due  to  the  wing 

and  include  the  perturbing  effect  of  the  free  stream,  we  can  write 
Equations  (24)  through  (26)  as 


Equations  (27)  through  (29)  define  the  coordinates  of  any  point  in  the 
wake  mesh  under  the  influence  of  wake  contraction  (axial  and  radial 
velocity  variations),  wing,  nacelle,  and  free-stream  velocity. 

In  the  above  development,  it  has  been  assumed  that  the  blade 
azimuthal  increments,  £<//  ,  are  all  equal.  The  radial  increments  may 
have  any  distribution  desired. 


<- 
Wx  _ 

?sj  - 

_  WY  _ 
" 

Wz  _ 

?SJ  ~ 
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LIST  OF  SYMBOLS  FOR  APPENDIX  I 


% 

I* 

% 


*>X 

m 

V 

*r 


£  integer  denoting  azimuthal  position  in  the  wake 

J  integer  denoting  radial  position  in  the  wake 

k,  integer  denoting  radial  position  of  a  blade  element 

KA  integer  denoting  azimuthal  position  of  blade 

a 7a  total  number  of  azimuthal  stations 

A4?  total  number  of  radial  stations 

X  component  of  velocity  at  jth  wake  point  at  time  S 
due  to  the  wing,  ft/ sec 

Y  component  of  velocity  at  wake  point  at  time  6 
due  to  the  wing,  ft/ sec 

Z  component  of  velocity  at  J**  wake  point  at  time  S 
due  to  the  wing,  ft/sec 

radial  distance  from  centerline  of  rotation  to  the  jih  wake  point,  ft 

distance  from  centerline  of  rotation  to  the  blade  tip,  ft 

free-stream  velocity,  ft/sec 

radial  velocity  of  the  ring  at  time  6,  ft/sec 

radial  velocity  at  jth  wake  point  at  time  6 
due  to  the  nacelle,  ft/sec 

number  of  revolutions  of  wake 

axial  velocity  of  the  ring  at  time  5,  ft/sec 

axial  velocity  at  jih  wake  point  at  time  S 
due  to  nacelle,  ft/ sec 

nonrotating  coordinates  fixed  at  prop  hub 
coordinates  of  the  ij*h wake  point 
Z  angle  free  stream  makes  with  thrust  axis,  deg. 
f  azimuthal  position,  deg. 

Af  azimuthal  increments,  deg. 

-CL  rotor  speed,  radians/sec 


n 


rj 


% 

Vl, 


JY 

*6t 


<X..Y;Z*) 
(£ >  £kj 
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APPENDIX  II 


THE  WING  AND  NACELLE  MODELS 


NACELLE  MODEL 


To  model  the  aerodynamic  effect  of  the  nacelle,  we  replace  it 
by  a  Rankine  body  of  equivalent  length  and  projected  frontal  area. 

Since  the  velocity  distribution  induced  by  the  Rankine  body  is  axially 
symmetric,  only  the  radial  and  axial  velocity  components  need  be 
given: 


VN  (rtl) 


h  -{h*+a,1 

X  X 

(30) 


WH(r.z)* 


Vffh-J h  *+(Ll 


(L  -  a,  +  z) 

U  +  cL  +2) 

y.U-ei  +  Z)‘]* 

[r*f  (L  +di+Z)2Y/z 

(31) 


where 
VN  = 
W*  - 

K?  = 

h  = 
L  = 
a>  = 

X  = 


radial  velocity  induced  by  nacelle 
axial  velocity  induced  by  nacelle 
velocity  which  the  nacelle  "sees" 
Rankine  body  radius  (maximum) 
Rankine  body  semilength 
~  (note:  to  determine  5,, 

L_ 

h 


must  be  solved ) 


P  =  radial  position  at  which  velocity  is  desired 
Z  =  axial  position  at  which  velocity  is  desired 
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If  we  consider  the  velocities  yvand  W^to  lie  in  a  plane  defined  by 
the  axis  of  rotation  of  the  prop  and  the  blade  radius  at  some  azimuthal 
position,  f  ,  then  we  can  resolve  the  l^and  Wvvelocity  components 
into  the  fixed  (Xa,  Y0  ,  Zc  )  system  shown  in  Figure  19  .  We  need  only 
multiply  V*by  cosy  and  sm^to  obtain  the  velocity  components  along 
X0  and  Y0  axes,  respectively.  The  (^component  lies  along  the  ZQ 
axis  and  needs  no  further  manipulation. 

Note  that  the  (  V*,  wM)  of  this  Appendix  correspond  directly  to 
the  (  Vg\  of  Appendix  I.  The  subscripts  Sj  denote  the  evaluation 
of  (  at  the  points  in  the  wake  identified  by  S  and  j  . 

To  obtain  the  preferred  to  in  Equation  (8),  page  12,  the  ( 
are  evaluated  at  the  jth  blade  element  and  resolved  perpendicular  to  the 
airfoil  chord. 


WING  MODEL 

The  wing  is  represented  aerodynamically  as  a  bound  vortex 
placed  at  the  quarterchord  of  the  wing.  The  velocities  induced  by 
the  wing  at  any  point  in  the  field  are  given,  in  the  fixed  coordinate 
system,  as 

<  [m»”b  (32) 

$"r  *  <33> 

<r  -  £  *  [4-b  -  "vvi  <3-» 
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where 


r#  -  wing  bound  vortex  strength 


K  = 


rarb[ra.rb 

l  s2  2  2  \ 

ri  *  Ji  +mi  +  ni 


<-  x-*i 

mim  y-  m 
*-G- 


f  i  •  a.,  b 


(X  ,Y  ,Z  )  coordinates  in  the  fixed  system  of  the  point  at  which 
the  wing  induced  velocity  is  desired 

a"  end  of  bound  vortex 

"b"  end  of  bound  vortex. 


The  velocities  (  y*" ,  ,  y*z)  correspond  directly  to  the 

(  fa*  »  (ffl  of  Appendix  I.  The  subscripts  <5j  denote  evaluation  of 
(y"">  f*ri  y.w* )  at  the  points  in  the  wake  identified  by  <$  and  j  . 

Similarly  the  of  Equation  (8),  page  12,  is  obtained  by 
evaluating  the  velocity  components  of  this  Appendix  at  the  jth  blade 
element  for  a  given  blade  azimuthal  position  and  resolving  them  into  a 
component  normal  to  the  airfoil  chord. 

See  Appendix  III  for  more  information  on  the  derivation  of 
Equations  (32),  (33),  and  (34). 
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LIST  OF  SYMBOLS  FOR  APPENDIX  II 


a, 

h 

K 


x 

z 


(note:  to  determine  a,  p*  ^  z-jf+a*  must  be  solved) 
Rankine  body  radius  (maximum),  ft 

[r<t*  rb\  /  ra.  rb  [r*  rb  *-4  -*6  +  mOLmb  +  nd  "b] 

Rankine  body  semilength,  ft 

*-?*>. ft 
•  ft 

,  ft 

X  component  of  the  wing  induced  velocity,  ft/sec 
Y  component  of  the  wing  induced  velocity,  ft/sec 
2  component  of  the  wing  induced  velocity,  ft/ sec 
radial  position  at  which  velocity  is  desired,  ft 

*  ”1  »  ft;2 
2 

,  ft 

radial  velocity  induced  by  nacelle,  ft/ sec 

velocity  which  the  nacelle  "sees”,  ft/ sec 

axial  velocity  induced  by  nacelle,  ft/ sec 

coordinates  in  the  fixed  system  of  the  point  at  which  the 
wing  induced  velocity  is  desired 

L_ 

h 

axial  position  at  which  velocity  is  desired,  ft 

2 

wing  bound  vortex  strength,  ft  /sec 


(**•*?*'£*)  "a"  en<*  ^ound  vortex 
( %b>  *}b>  "b"  end  of  bound  vortex 
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APPENDIX  III 


DERIVATION  OF  VELOCITY  INDUCED  BY  THE  WAKE 

The  induced  velocity,  £  ,  at  any  point  P  due  to  an  arbitrarily 
oriented  vortex  element  of  strength  /  per  unit  length  is  (in  vector 
notation)  v  ^ 

**  \w*iWwfr)  <35> 

n  7  fbdsxr 

?  =  V/r^  If  I*  (36) 

where 

<^5  -  *  +  cLYj  +  atZ^h 

T, J.k  =  are  unit  vectors  along  the  axes 


K 


SENSE  OF  /  POSITIVE  FROM  a  TO  b. 


Figure  20.  ARBITRARILY  ORIENTED 
VORTEX  ELEMENT 

The  velocity  components  of  Equation  (36)  can  be  written  in  terms  of 
the  coordinates  of  the  end  points  of  the  elemental  vortex  and  the 
coordinates  of  the  point  at  which  the  velocities  are  desired.  Thus, 


r  +  ir  1 

ft!  \  n  ^  ^ 

+TT  L rarb  (rar*  *  n anb)  J[  J 


(37) 
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f  '*«•. 

<V 


J 


£**  ’  +Tr\reLrb(r<Lrb  -hUCb2b  +  m^b  +  /ia.nfr][n«''^  "^wt>] 

qM  =^r - - - ^4^ - Ik^-Ul 

*  °  ^  L  r*rbf r4  rb  *Jaib*  mamb  ♦  na  nb  J  |_  a  J 


(38) 


(39) 


where 

2  2  2 

ra,  *  A^mcx  +'”a. 

t  .  I  |  2 

rb  ~  '“b  *  mb  *  nb 

-4.*  *-*•. 

^  ^ 

*Y  *  y/-7a, 

*  y'-/^ 

'V’z'-r«, 

^  •  *#-a . 


The  components  given  by  Equations  (37)  through  (39)  are  applicable  to  both 
the  bound  vortex  of  the  wing  and  the  elemental  vortices  of  the  wake. 

To  determine  the  induced  velocity  perpendicular  to  the  blade 
chord  at  any  azimuthal  and  radial  position  of  the  blade  due  to  a  single 
arbitrarily  oriented  vortex  element,  we  must  make  a  transformation 
from  the  fixed  (Xe  ,  Y0  ,  lD)  coordinate  system  to  the  (X,  y,  z  )  coordinate 
system  rotating  with  the  blade.  Consider  Figure  21  (all  quantities 
shown  are  positive). 
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Figure  21.  INDUCED  VELOCITIES  IN  FIXED  AND  ROTATING  FRAMES 

Thus,  in  terms  of  the  components  )  and  the  azimuthal 

position  of  the  blade,  the  perpendicular  velocity  can  be  written  as 


s  +y.t'Cose. 


If  we  now  define  the  geometric  portion  of  Equations  (37)  through  (39)  as 
follows, 


(ra,  +  r>) 


Y  U  WtCra,  r6  +  m*rrn+  n^)  (' )  "  °it 

<}Yo  _±_  \ _ ^  +  rb) _ I 

Y  ^  /a rb(ra.rb-h^<L^k+rnournb-hr,a.^b)  _  ^ 


I**  _l  r 

r  ~  VW  ra-rb  (rCLrb 
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then  Equation  (40)  can  be  written  as 


^  '  [{ °le  ^  "  °Y0  £°s  &  }  5cn  &  *■  cos  y 

or  simply  as 

'U'1  -  5^/. 


(41) 


Equation  (41)  was  written  for  one  point  and  one  vortex  element.  To 
determine  the  net  contribution  at  some  point  j  due  to  all  the  vortices 
i  in  the  wake,  we  simply  form 


O 0 

L 

is  f 


(42) 


where  denotes  the  geometric  influence  function  of  the  ith  vortex 

element  at  the  jth  point. 
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LIST  OF  SYMBOLS  FOR  APPENDIX  III 

^  ^  ^  » 

unit  vectors  along  (X0.YO,Z0)  axes 

4 .*  *'-!*.»  ft‘ 

ma.b  Y'-fa.f  ft 

nA.  b  *  *  %a.  6  •  ^ 

(fxtt 9*1*9*  )  wa^e  induced  velocity  components  referenced  to  {X0,Y0tZ0) 
system,  ft/ sec 

*&,b  *  ml  b  +  n+  b  »  ft2 

t*. 

V'j  velocity  perpendicular  to  the  airfoil  chord  of  the  j  blade 
element  due  to  the  wake,  ft/sec 

(X,y,z)  coordinate  system  rotating  with  the  blade 

(X,  Y',Z')  coordinates  in  (Xs,  Ya  ,  Za  )  system  of  point  at  which  wake 
induced  velocity  is  computed 

(Xs,  Y0,Zt)  nonrotating  coordinates  fixed  at  hub 

"Yi  strength  per  unit  length  of  the  vortex  element,  ft  /sec 

6  angle  between  the  airfoil  chord  and  the  prop  plane,  deg. 

(<XZU,  crr .  induced  velocities  in  ( X .,  Y,  ,  Z0  )  directions  due  to  a  vortex 

element  of  unit  strength,  1/ft 

geometric  portion  of  the  Biot-Savart  law  for  the  effect  of 
the  ith  vortex  element  at  the  jth  point,  1/ft 

f  azimuthal  position  of  the  blade,  deg. 
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APPENDIX  IV 

DERIVATION  OF  THE  ELASTIC  EQUATIONS  OF  MOTION 

A  Lagrangian  approach  is  used  to  derive  the  equations  of  motion 
of  the  rotating  propeller.  In  addition  to  the  previously  defined  coordinate 
systems,  (X0  ,  Y0  >  70  )  fixed  in  space  and  (  X  ,  Y  ,  Z)  rotating  with  the 
blade,  the  local  coordinate  system  (  % ,  rj  ,  $  )  as  shown  in  Figure  22 

is  introduced.  In  Figure  22  the  axes  5  and  X  are  normal  to  the  plane 
of  the  paper. 


Figure  22.  LOCAL  BLADE  SEGMENT  COORDINATE  SYSTEM 


The  airfoil  is  assumed  to  be  thin  and  symmetrical.  The  components  of 
the  total  velocity  in  the  rotating  frame  of  any  point  P  on  the  airfoil  is 
given  by  (dot  indicates  differentiation  with  respect  to  time) 

Vx  =  X  -  ClY  +  Vf  sin  T cos  (43) 
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VY  =  f  +  flX  -  Vf  sin  T  cos  4 


(44) 


Vz  s  Z  *  Vf  cos  V 


(45) 


where 


X  -  % 

Y  »  eQ  +  *]  cos  6  ~  £sLn  O  *  trcos/O  -  u/sinfl 
Z  =  rj  sin  9+  gcas  9  *  urcos^Q  +  U'Scn  jQ  , 


and  the  displacements  (v-yur)  are  taken  parallel  and  perpendicular 
to  the  undisturbed  airfoil  chord.  The  displacements  )  are 

expanded  in  terms  of  the  generalized  coordinates  (  ur(P\  irW,  4>(m)  ), 
which  are  the  tip  deflections  of  the  uncoupled  flapwise,  chordwise,  and 
torsional  vibration  modes  of  the  blade: 


ur*  y 

(46) 

t 

u- £  v'Vf'V 

(47) 

*>  ■■  z*™*? 

(48) 

s 


where 

-  flatwise  mode  shape  of  the  t mode 


70 


fj?)s  edgewise  mode  shape  of  the  mode 
torsional  mode  shape  of  the  s*h  mode 

tip  deflection  of  the  flatwise  mode;  dependent  on  time  only 
tr<P- tip  deflection  of  the  ^  edgewise  mode;  dependent  on  time  only 
<p(s)-  tip  deflection\>f  the  sth  torsional  mode;  dependent  on  time  only. 


At  this  point  it  was  decided  to  employ  the  coupled  flatwise-edgewise 
modes  as  the  generalized  coordinates.  Hence,  Equation  (47)  became 


(49) 


ir  =  ^_  ur' 

P 

where  -  constant  determined  from  a  coupled  mode  analysis  as  is  the 
coupled  mode  shape  4^  .  Equation  (46)  has  the  same  form  when  the 
coupled  mode  shape  4^  is  used. 

Lagrange's  equations  of  motion  are  given  as 


oL  /dT  \  dT  _  qL  fdU  \ 

d?r  \a?r) 


du  so  „ 

4  “z -  +  \ ",  ■  s  Qr 

<*?r  ' 


(50) 


where 


T  =  kinetic  energy  of  the  system 
U-  potential  energy  of  the  system 
D-  dissipation  function 
Qr  -  generalized  force 

generalized  coordinate 
=  r*h  generalized  velocity 

r  =  number  of  generalized  coordinates  required  to  define  the  system. 

The  hub  is  assumed  to  have  infinite  impedance,  and  therefore  each  blade 
may  be  treated  as  an  independent  cantilever  blade.  Hence,  the  total 
kinetic  energy  of  one  blade  is  given  by 

T~if  (v‘*v**vz)dm  <51> 

2  vauet 

M MSS 
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where  for  a  blade  of  uniform  density, 


dm  -  pdAdX. 


The  potential  energy  is  (following  Reference  9) 

u- if  ]$)-*£**££!** 

r  u 


where 


* 

X  A®  dX 


dfi_ 

dX  dLX 


t 

A 


cL*ur  dp  dtr  d*P 

dx*  +  2  dx  ~dx  +  d}Tl  ^  '  \  ax) 

„2cLfi_^_c£L a, -(&)*„ 

dX*  dX  dX  dX  \dX J 


d<P  cLus  d  *4>  dfi  oL4> 

dX  dX  ’  dX  *  ^  ‘  d  dX  dX 


X 


rr>X  dX 


(52) 


£If  =  flexural  rigidity  about  major  neutral  axis 
Eli-  flexural  rigidity  about  minor  neutral  axis 
GJ  =  torsional  rigidity 

e»  *  ifA  idA 

k‘  -  jJa  n*A 

EB ,  «  E^  q2(i2-k,j)  dA 
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£Bt  -  ej nC^1  -kl)d.A. 


The  dissipation  energy  is  obtained  by  assuming  that  the  total  system 
damping  can  be  defined  in  terms  of  a  summation  of  the  damping 
proportional  to  the  displacement  and  in  phase  with  the  velocity  of 
each  degree  of  freedom  of  the  system  when  they  are  individually 
performing  simple  harmonic  motions. 


ase 


The  tip  deflections  are  expanded  in  a  Fourier  series  on  p  , 
azimuthal  position  (which  is  equivalent  to  time  for  the  steady-state 
e):  _  _ 

1L'  f;' 

=•  a/(J£  +  2^  cos(np)+2  b(J?n  sin(nP) > 

nwt  n*i 

<p(s)  *  <2^  +  /  cos(np)  +  £  bp*sin(np). 

n*i 

generalized  force  is  similarly  expanded  in  a  Fourier  series: 

Ei-t 

y  c.os(n<£)  +■  /  sm  (n <k) . 

n:/  n*f 

f-t  < 

-  <>+Z  " 


(53) 


(54) 


(55) 


The  Lagrangian  is  applied  to  Equations  (51)  through  (56), 

and  like  coefficients  are  equated.  The  resultant  set  of  equations  is  given 

in  the  following:  For  the  0*^  harmonic,  n  -  0 


1 

f 

'i 

-n2 

I 

r\l  pU’j) L 

11  [L*’t 

„(j) 

CO<*o 

1 

1 

< 

,  -  < 

- n 2 

1 

l\ 

n‘  [e?j>]  <  rj‘J> 

4>o 

_  g-(i>  _  n*A(*' 

1 

J 

^  > 
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For  1st  to  harmonic ,  n  =  1,  2,  ... 


{NA/2  -  / ) 


-nWeUP+F™  j 

\  &  Fj?> 

1 

<j)  r(ij) 

~St  G^t 

(a.1*' 

n 

1 

1 

1 

1 

0  /I  ^0ur  ~  ^4ur  1 

I 

V/l*sJjV  F(ji} 

c?i> 

/>  *j; 

„  U) 

*+n 

K‘ 

-A*  ]C'&  -  P%> j 

1 

1 

1 

1 

1 

1 

►  *< 

► 

1 

_  O)  r?  <*J)  | 

9uJ-  ’ur 

r  <*j) 

1 

1  -ntl'a™*  F^> 

,  (J) 

b^n 

Vi) 

l{4^» 

n 

\ 

1 

1 

1 

1 

1 

(J)  r(ij)  | 

tyov'  tM' 

rr  <  ij) 

9+  FP 

\-«WD$H  -G%> 

-nri’df'+F^ 

,<J) 

bt 

Tn 

fc‘> 

*V> 

\ 

1 

1 

1 

1 

+  n*  Ep^ 

- 

/ 

where 


S  f  r»(eo  + C  Cos Stnfl  -  cos/i)  d.X 

ro 

P  R 

A*f  =  /  (me0e +p.  cos/6)  sin J& -f^dX 
°r0 

c* =  cvrcvr;  ^ 

#*•  /Vv;'v, 

"o 

Cl^p  J  (me  sm /3)(fJ?} sinji  -  cos fl)  -F^ oiX 

ro 

D!%  s  J'me&W'dX 

Ur, 

=  J  (me0  t  cos  £  +  cos  2fi)  -F^-ffl^ctX 

ro 

F(Sp)=  J  (ei,  tf'-Fr”  +  EIZ  ^r>fip))  OCX 

9  o 

Fj?  *  Tie,  f?ctX 

re 

F“>-  f*  T  l*  d?  Ct**‘  M 

F' *  =  J.  z  k‘  dT  Sx~  d* 

9  0 
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o’  -  ffafrr  f  -  z"*£?)** 


Q<re) 

^  us 


s  J''  mx9 


(rp'i 

us 


cLX 


r~(e 


rtf, f  .»*”)  tjfM 


i  r* 

X  S  INITIAL  AXIAL  TENSION  =  fl  /  mX  dX 

*a 

_  cdVjJ1  ^  dfi  .(p)  df/  d  fi  (,)  (f)  /dfi\ 

ff  *  ~dx*~  * 2  ^7  v  ~M~  7x*  kv  fr  (277' 


fir)  = 


Tus+ 


=  6 


=  2 


^V,‘rt  df^  d‘fi 

~dtl 


r(N 

'US 


cLX  dX 


3  r(f>  f^\  L(rt r(l 

UyJ 

dj£  dC' .  *  2  2^  km  r, ft 

J~  -•  ”  dx  dX  *  * 


dX  oLX 


Q<rp)  _ 
<3  or 


Jm  l \dx  dx  */\dx  ax 


(A) 
is  I 


ftu)  dff  dfi  Pm\(.m  df”  d£  -in )]*y 

(r~-3T  ‘ ■3TfWr*'  ~dT‘  dx^)\ax 
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LIST  OF  SYMBOLS  FOR  APPENDIX  IV 


n,f>)  A <f>) 

&  ,bu,n 

Fourier  coefficients  of  blade  flatwise  response  in  the  p**> 
degree  of  freedom 

A 

blade  cross  sectional  area 

n(S)  bf9) 

th 

Fourier  coefficients  of  blade  torsional  response  in  the  i 
degree  of  freedom 

D 

dissipation  function 

£3, 

Ef*r)*(n2-kt)dA  ,  ft4 -lb 

E3Z 

Ej[n(nt-^)dA  .  ft3-ib 

EZ, 

2 

flexural  rigidity  about  major  neutral  axis,  ft  -lb 

EIZ 

flexural  rigidity  about  minor  neutral  axis,  ft^-lb 

e 

distance  from  elastic  axis  to  center  of  gravity  of  section,  ft 
((e=  0)  for  the  case  shown  in  Figure  22) 

ndA  > £t 

*-o 

distance  from  the  K  -axis  to  the  elastic  axis  of  the  undisturbed 
airfoil,  ft 

c  <p)  S-  (p) 

Fourier  coefficients  of  the  generalized  force  in  the  pth 
degree  of  freedom 

%<*) 

*  <^n 

Fourier  coefficients  of  the  generalized  force  in  the 
degree  of  freedom 

f(V 

'  *r 

edgewise  mode  shape  of  the  ^^m  j 

r  0*> 

’  ur 

flatwise  mode  shape  of  the  pth  mode 

■fiS) 

torsional  mode  shape  of  the  mode 

GJ 

2 

torsional  rigidity,  ft  -lb 

ki 

if/dA,  ft2 

LA 

flatwise  to  edgewise  coupling  constant  in  the/^mode 

m 

blade  mass  per  unit  span,  slugs/ft 

Qr 

generalized  force 

9r 

generalized  coordinate 

?r 

/•^generalized  velocity 
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R  radial  distance  from  centerline  of  rotation  to  blade 
tip,  ft 

r0  radial  distance  from  centerline  of  rotation  to  blade 
root,  ft 

T  kinetic  energy  of  the  system 
u  potential  energy  of  the  system 
free-stream  velocity,  ft/ sec 

v ’  elastic  deformation  of  the  elastic  axis  parallel  to  the  undisturbed 

position  of  the  airfoil  chord,  ft/ sec 

tip  deflection  of  the  qt*  edgewise  mode;  dependent  on  time  only,  ft 

Vx>  VY  ,  Vz  ( X,Y,Z  )  components  of  the  total  velocity  of  any  point  on  the 
blade,  ft/ sec 

or  elastic  deformation  of  the  elastic  axis  perpendicular  to  the 
undisturbed  position  of  the  airfoil  chord,  ft/sec 

tip  deflection  of  the  ^  flatwise  mode;  dependent  on  time  only,  ft 

(x.Y.Z  )  rotating  coordinates  fixed  at  the  hub 

(Xe,Ye,  2.)  nonrotating  coordinates  fixed  at  the  hub 


/i  angle  of  built-in  twist  with  respect  to  the  prop  plane,  radians 

0  total  geometric  angle  of  blade  section,  radians 

jlu  section  polar  mass  moment  or  inertia,  slug/ft 

({•If  £ )  local  coordinates  fixed  at  the  blade  elastic  axis 

3 

/ o  blade  local  density,  slugs /ft 

7  angle  free-stream  velocity  makes  with  thrust  axis,  deg. 

<p  torsional  displacement  of  a  given  blade  element,  radians 

<flt)  tip  deflection  of  the  s*A  torsional  mode;  dependent  on  time  only, 
radians 

V  azimuthal  position  of  blade,  deg. 

71  rotor  speed,  radians/ sec 
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APPENDIX  V 


DERIVATION  OF  THE  GENERALIZED  FORCES 


By  definition, the  product  of  the  generalized  force  timea  the 
generalized  coordinates,  summed  overall  degrees  of  freedom,  must 
1  the  total  virtual  work  done : 


II  Fi  *  d ri  =  X 


where 

-  net  force  acting  on  particle  (vector) 

Qj  =  jth  generalized  force 
oLr c  -  displacement  oi  ith  particle  (vector) 

=  differential  of  j**  generalized  coordinate 
i  -  number  of  particles 
c J  -  number  of  generalized  coordinates. 

Now  since  the  position  vector  r-  is  a  function  of  the  Qj's  , 


hence, 

For  each  j  ,  therefore, 


Zr, -(§)••>. 


(57) 


Now,  to  form  the  virtual  work  done  by  the  force  acting  on  the  prop 
blades  due  to  a  virtual  infinitesimal  displacement  of  the  generalized 
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coordinates,  define  the  following  geometry: 


Figure  23.  GEOMETRY  FOR  DERIVATION  OF  GENERALIZED  FORCES 


where 

lift/unit  span  =  /,  taken  positive  perpendicular  up  from 
local  relative  velocity 

drag/unit  span  =  cf,  taken  positive  in  direction  of  local 
relative  velocity 

moment/unit  span  =  m ,  taken  positive  leading  edge  up 
about  the  elastic  axis. 

Then  at  each  section  of  the  blade  and  each  azimuthal  station, 

J  s  ^0  V r  -  bp  V*Cj  (58) 


d  =  bp  V*cd 


(59) 
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m  -  e  [j cos  at  *  dLsin  a,]  *  2b‘fi  V *em 


(60) 


Now  the  virtual  work  is  given  by 

R 


dW  «  J  [fxoLX  +  FYdY  +  FzdZ  +  moie}cLX 


(61) 


where 


r-  c  dur, 

F*  mju*  -dx^-dT 


Fy  ~  CO «  j  sir)  °*X  ~  [rf]  coS  •=»”^  .»*/■>  Qfx  •+  ct  COS  &x  J 

Ft  =  | yZcos  ]  co5  arx  ~  \^icos  ccx  ~dsm  ax  j 


<xx  -  arc  tan 


:$] 


Vz  -  total  velocity  perpendicular  to  (X,Y  )  plane 
Vy  -  total  velocity  in  the  plane  of  [X  ,Y  )• 

To  obtain  the  differentials  dX .  e/Y,  dZ,  a.ndd&  e  the  following 


ig.e*i„9-3-2eco,o 

dX  dX 


(62) 


*  ‘  e0  +  v; 


(63) 


Z-  ur4 


(64) 


81 


where 


v',  -  ircosfi  -  cv'scnfl 
u^t  -  ur  cos  fi  +  xr  sin  j5 


The  displacements  {  u- ,  ur ,  <p  )  are  expanded  in  terms  of  the  generalized 
coordinates  (a^p\  u-(Pl  <p(f))  which  are  the  tip  deflections  of  the  flapwise, 
edgewise,  and  torsional  vibration  modes  of  the  blade: 

P 

ir  =  ^z^(p)4p) 

p  p 

<P  -  Z  • 

* 


Introduce  the  small  angle  approximations 

(p  «  / 


and 


sin  QJe  ~  oce 


Form  the  differentials  of  Equations  (62)  through  (64)  and  substitute  into 
Equation  (6 1)  to  obtain  the  virtual  work: 

-  Z  J  (Jan  ax  +  <£cos  az  )(kl^  f^cos/B  -  dXdur  (p) 


?  r> 


Z  J  (dcos  ax  -  &  sin  ax )  (fZc  oS  +  *'£-f")sinfi)dXdur‘p) 

?  ro 

Li  W  cos  <xe  +<£sih  cce)  +  Zp  V  zb2cm  -F^dXd^V 

Q 


(65) 


Apply  Equation  (57)  through  Equation  (65)  to  obtain  the  generalized  force: 


.  (ft  f  (r)  dfi 


us  ,  l" yr<  — 

"  oivj 


“W**- -^ls‘"4dx 


rR 

-J  (jsir)  ax  +  d  cos  az)( c<?s  yS  -  +£?.sLn  fi)  d* 


(66) 


and 


/  (V  c<75  <ZX  -d  sin  a  j)  (■?"*  cos/3  +  k^-P^sinfit)  cLX  r  *  ^  2 ^ 
r0 


(5^  ^  =  J*  jgf^cos  +  d  Sen  Qte}  4- 2^>V  (?  dX%  f  -  1,2.  (67) 


Note  that  Equations  (66)  and  (67)  correspond  to  the  left-hand  side  of 
Equations  (55)  and  (56)  of  Appendix  IV,  page  73. 
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LIST  OF  SYMBOLS  FOR  APPENDIX  V 


b 

Cat 

CJ 

C”c/¥ 

cL 

d9j 

Co 

t(r) 

Tir 

L  <P) 
rur 


Jt 

m 

Qa 

Y 

J'V-  Vz 

'if,  US' 


If,.  US', 


uy(P) 

(X,Y,Z) 

<*e 

az 

A 


local  blade  semichord,  ft 
section  drag  coefficient 
section  lift  coefficient 

section  moment  coefficient  about  quarterchord 
drag  per  unit  span,  lb/ft 
displacement  of  particle  (vector),  ft 
differential  of  j ^generalized  coordinate 

distance  from  X  -axis  to  the  elastic  axis  of  the  undisturbed 
airfoil,  ft 

net  force  acting  on  i**  particle  (vector),  lb 
edgewise  mode  shape  of  the  pih  mode 
flatwise  mode  shape  of  the  pn  mode 
torsional  mode  shape  of  the  mode 
number  of  particles 
number  of  generalized  coordinates 

flatwise  to  edgewise  coupling  constant  in  the  pth  mode 

lift  per  unit  span,  lb/ft 

moment  per  unit  span,  ft-lb/ft 

^generalized  force 

total  local  velocity,  ft/ sec 

components  of  the  total  velocity  in  the  plane  of  (X,Y)  and 
perpendicular  to  the  plane  of  (X,  Y),  ft/sec 

elastic  displacement  of  the  elastic  axis  parallel  and  perpen¬ 
dicular  to  the  undisturbed  airfoil  chord,  ft 

elastic  displacement  of  the  elastic  axis  parallel  and  perpen¬ 
dicular  to  the  (X,  Y  )  plane,  ft 

tip  deflection  of  the  p*h  flatwise  mode,  ft 

rotating  coordinates  fixed  at  the  hub 

local  effective  angle  of  attack,  radians 

angle  local  total  velocity  makes  with  the  prop  plane,  radians 
local  built-in  twist,  radians 
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e  distance  from  elastic  axis  to  quarterchord,  ft 

O  total  geometric  angle  of  blade  section,  radians 

3 

/o  density  of  air,  slugs/ft 

<p  total  elastic  torsional  displacement  of  the  airfoil  section, 

<P(S)  tip  deflection  of  the  torsional  mode,  radians 


radians 
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APPENDIX  VI 


DERIVATION  OF  BLADE  STRESSES 


From  Reference  7,  the  longitudinal  strains  were  derived, 
Assume  that  the  airfoil  consists  of  a  bundle  of  fibers.  Now  pass 
a  plane  perpendicular  to  the  elastic  axis  (EA)  and  consider  the  location 
of  a  fiber  through  the  point  P  before  and  after  a  deformation: 


Figure  24.  ELASTIC  DEFORMATIONS 


86 


X  +  - 


aittf 

dX 


cos  9  -  £sin  ^qslnd  +  Ceos  (68 ) 


Y  +  4-  ["17  COS  G  -  Csin  @]  +  [17 


C/M  ^  -A.  & 


(69) 


Z,  -  Z 

C  ui  >vi‘  “}9  = 

0  = 
(x.y.z)  = 
(x,X.A)  = 


+  u^  +  y\sinG  *  Ceos  0J  . 

displacement  in  (X  V,  Z) directions  of  point  P  from 
its  undeformed  position 

built-in  twist  +  torsional  deformation  =  A*~<p 
coordinates  of  the  point  Pt  prior  to  deformation 
coordinates  of  the  point/3,  after  its  deformation 


(70) 


To  find  the  change  in  fiber  length  of  the  fiber  through  the  point  P  ,  form 

(ok,)*  *  (dX,)*  +  (dr,)*+  ( dZ y.  (71) 

The  original  fiber  length  was 

( ds )*  =  (dX)*  +  (dY)z  +(ail)2'  (72) 

The  tensile  strain  in  the  fiber  is  defined  as 

e  )  ~  (ds)  t  (73) 

ds 

Thus  (the  primes  indicate  differentiation  with  respect  to  x  ), 


€  -  u.'f~  r\(y?cosfl  +  W,'sin/3  )  -  C(~  Vfstnfl*  W,"cos jS )  *  (r\3+  C  *) (74) 


or  in  terms  of  the  displacements  parallel  and  perpendicular  to  the 
undisturbed  airfoil  chord  (v-,  1^),  introduce  the  transformation 


87 


vJ1  -  trcos JS  -ursCnfi 


(75) 


■  t/'Scn  J3  -  or  cos  jt  . 


(76) 


Since  the  fiber  is  assumed  to  be  effective  in  uniaxial  strain  only,  the 
stress  follows  immediately  as 

o~  =  Ee.  (77) 


To  eliminate  uj,  equilibrium  of  forces  in  the  X  direction  is  imposed: 

X  -  E  I  edfd/i  (78) 

VOrt  -t/2 

where 

( Injure)-  leading-  and  trailing-edge  coordinates 
t-  the  airfoil  thickness 


X  =  centrifugal  and  aerodynamic  forces  acting  in  X  direction 


R 

X  fL *  dLm 


The  assumption  of  a  thin  symmetrical  airfoil  is  applied  and  the  equation 
is  solved  for  aj ;  the  result  is  used  in  Equation  (74).  Now  the  transformation 
of  coordinates  is  made  and  Equation  (77)  becomes,  in  terms  of  the 
displacements  parallel  and  perpendicular  to  the  chord, 


-  r  [w*+  tjs'v'-t p'v -(j5')*wJ 


(79) 


where  the  primes  indicate  differentiation  with  respect  to  X  and 
E  -  Young's  modulus 
A  =  airfoil  cross-sectional  area 


(for  a  symmetric  airfoil). 


To  account  for  the  torsional  strains  and  stresses,  it  is  assumed  that 
there  exists  a  St,  Venant's  torsional  resisting  stress  which  would 
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develop  if  the  blade  were  not  initially  twisted,  plus  a  shearing  component 
developed  by  the  longitudinal  stress  in  the  fibers.  This  component 
develops  because  of  the  live  twist,  ,  and  also  the  built-in  twist,  /5  . 


Figure  25.  PICTORIAL  REPRESENTATION  OF  TORSIONAL  DEFORMATIONS 

The  in-plane  component  of  fiber  stress,  a -,  is  (see  Figure  Z 5) 

crs  =  crsin  y  s;  a~  y, 


but 


thus, 


c  •  -  ffj  ■ 


The  total  resisting  torque  is  given  by 


r  iLf  r  i/2  - - 

q  =  gjq  +  J  J  -| in *+ 


(80) 


he  p  i/2 
Hre  -t/z 

where,  after  accounting  for  a  thin,  symmetric  airfoil,  Equation  (81) 


(81) 
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becomes 


« -[w  «*££)' 


(82) 


where 


0  -  shear  modulus 


B,  *  f  rj2(n*'  k%)  dA 
JA 

&z-  JAr[(r\* '  ^ 

jJa1Ua 


A  -  -  2p’ur'~ &"ov  + 

ZLy  =  area  moment  of  inertia  about  £ 
A  =  cross-sectional  area 
c  =  sectional  chord 


To  determine  the  maximum  shear  stress, the  section  was  treated  as 
though  it  were  in  pure  torsion.  Then  for  a  symmetrical  airfoil  (see 
Reference  1  2) , 


where 


°~s_ 


Qs 

J 


3 


D 


7T*  D 

16  A 


1  *  0.15 


7T ‘D* 
16  A* 


(83) 
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r  =  radius  of  curvature  of  boundary  at  point  of  tangency 
O  -  diameter  of  maximum  inscribed  circle. 

Thus,  in  terms  of  Equation  (82),  Equation  (83)  becomes 


majt 


*  7 [1&J +**&/}%  ^  + 


The  tensile  stress  at  any  point  in  the  cross  section  is  given  by 
Equation  (79), and  the  maximum  shear  stress  at  any  blade  section  is 
given  by  Equation  (84). 


(84) 
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LIST  OF  SYMBOLS  FOR  APPENDIX  VI 


A  airfoil  cross  sectional  area,  ft 

B,  J  r,‘tf-kf)dA  ,  ft6 

B2  Jn  !*!2-kS)dA  ,  ft5 
c  sectional  chord,  ft 

D  diameter  of  maximum  inscribed  circle,  ft 

/  2 

£  Young's  modulus,  lb/ft 


/  2 

&  shear  modulus,  lb/ft 

XY  area  moment  of  inertia  about  £ ,  ft 


*A 

Q 

r 


Ui'dA'*z 

total  resisting  torque  at  a  blade  section,  lb-ft/ft 
radius  of  curvature  of  boundary  at  point  of  tangency,  ft 


D  T  „  frr * D a  D  \]  ^ 

uTTlDl  L  ■  ('t A! 


1’ 

IT 

uS 

X 

(X,Y,Z) 

A 

r 


16A 

airfoil  thickness  distribution,  ft 

displacement  of  elastic  axis  parallel  to  the  airfoil  chord,  ft 

displacement  of  elastic  axis  perpendicular  to  the  airfoil  chord, 

centrifugal  and  aerodynamic  forces  acting  in  /  -direction,  lb 

rotating  coordinates  fixed  at  prop  hub 

local  built-in  twist,  radians 

angle  of  shear,  radians 

tensile  strain  of  a  fiber,  inches/inch 
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e 


local  total  geometric  angle  (=  radians 

J  V'*  -  2 ur  + 

C)  coordinates  of  a  point  P  in  the  coordinate  system 

fixed  at  the  blade  elastic  axis 

.  2 

cr-  tensile  stress  of  a  fiber,  lb/ft 

/  2 

<75  component  of  a~  in  the  plane  of  the  airfoil  section,  lb/ft 

(f>  total  torsional  displacement  of  a  section,  radians 


ft 
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APPENDIX  VII 


DERIVATION  OF  THE  QUASI -STEADY  DAMPING  TERMS 


To  obtain  the  quasi- steady  damping  terms,  follow  the 
development  given  in  Appendix  V,  page  79.  For  the  purpose  here,  define 
an  approximation  to  the  quasi-steady  lift  and  moment  at  the  blade 
section  as  follows:  (Dot  indicates  differentiation  with  respect  to  time.) 


4  " 


(QSD) 


(85) 

(86) 


Equations  (85)  and  (86)  are  used  to  form  the  virtual  work  from  which 
the  desired  quasi-steady  damping  terms  of  interest  are  obtained: 


Qj?>  .  +X Cr'*  «£'  *  Z  ^^7  (87) 

r  s 

-  -z -Z (88) 


where 


Q(J?  generalized  force  in  vr(?)th  mode 
generalized  force  in  <p^)th  mode 

RT  =  ffiba-C’l V.  co*  r(k(XHco*fi-C'^p) 
r.  1 

-  (m  XtZ'cosfi  *  d-x 

si%  *J pbo-n  v,c  f  V,  cos  r^'ty'cos p  -  p  ) 

-(Jfn  XC'cosjB  +  ks^-F^sin^dX 
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T‘Z'  «  H  cw  r)  *] \»t?dX 

(*0 

V<pfS)  3  J  />  bounce  €  |T(x/i)*+  (Vf  cos  r/j  f/vr^x  • 


The  above  quasi-steady  terms  are  added  to  both  the  right-hand  side 
and  the  left-hand  side  of  the  elastic  equations  of  motion  of  Appendix  III, 


(p) 


page  64,  and  treated  accordingly.  The  Q'"  and  will  appear  implicitly 


The  V^and  Vi***  win  aPPear 


in  the  J  ,  ,  and  Ja 

ft  /7  ' /? 

explicitly  as  elements  of  the  square  matrix  on  page  74. 


LIST  OF  SYMBOLS  FOR  APPENDIX  VII 


*/ 

*/ 

i  < p ) 

T  y. 

1  (p) 
'or 

s/t> 
k  p* 

y  (QSO) 

( QSD J 
m: 


UT<» 

X 

A 

cj 


*?Wj 


t 

P 

n 


section  lift  curve  slope  ( *<LcJi/dc[)4 

local  blade  semichord,  ft 

edgewise  mode  shape  of  p**  mode 

flatwise  mode  shape  of  mode 

torsional  mode  shape  of  ^  mode 

flatwise  to  edgewise  coupling  constant  in  mode 

approximation  to  the  quasi-steady  lift  at  the  blade  section,  lb 

approximation  to  the  quasi-steady  moment  at  the  blade  section, 
ft -lb 

free-stream  velocity,  ft/sec 

local  total  velocity,  ffVsec 

velocity  in  yp^bending  mode  at  blade  tip 

radial  distance  from  centerline  of  rotation,  ft 

local  built-in  twist,  radians 

distance  from  elastic  axis  to  quarterchord  at  blade 
section,  ft 

distance  from  elastic  axis  to  three -quarterchord  point  at 
J**  blade  element,  ft 

angle  free-stream  velocity  makes  with  thrust  axis,  deg. 
air  density,  slugs/ft^ 

velocity  in  <pih  torsion  mode  at  blade  tip 
rotor  speed,  radians/sec 
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APPENDIX  VIII 


EXPERIMENTAL  STRAIN  DATA 


The  purpose  of  mounting  strain  rosettes  on  the  blades  was  to 
obtain  the  true  state  of  stress  at  the  selected  blade  locations.  The 
strain  signals  from  each  gage  of  the  rosette  were  recorded.  Then,  at  a 
given  point  in  time,  the  three  strain  signals  from  a  given  rosette  were 
combined  to  obtain  the  principal  stresses  and  directions.  From  the 
principal  stresses,  the  longitudinal  and  shear  stresses  at  the  locations 
of  the  rosettes  were  determined.  The  reduction  of  the  measured  strains 
to  stresses  required  either  that  all  the  strains  be  read  simultaneously 
or  that  they  be  in  steady  state.  Neither  requirement  was  truly  satis¬ 
fied  for  these  data.  Some  of  the  strain  data  were  reduced,  however, 
by  reasoning  that,  if  any  record  could  be  found  for  which  the  propellers 
had  maintained  a  fixed  phase  with  respect  to  one  another  for  the  dura¬ 
tion  of  the  record,  the  strains  could  be  considered  periodic  for  that 
interval,  and  hence,  stresses  could  be  evaluated. 

The  mean  and  instantaneous  values  of  the  harmonically 
analyzed  experimental  stresses  and  strains  are  presented  in  Figures 
26  through  29. 

The  mean  longitudinal  strains  arc  otted  vs.  X ,  the  angle 
between  the  free  stream  and  the  thrust  axes,  in  Figure  30. 

The  effect  of  the  "clean"  vs.  "dirty"  configuration  on  the 
strains  is  indicated  in  Figures  28,  29,  and  30.  Figure  30  indicates 
that  the  "dirty"  configuration  had  caused  only  a  slight  increase  in  the 
outboard  (3.  67  ft)  mean  longitudinal  strain,  but  had  caused  a  three¬ 
fold  increase  in  the  mean  longitudinal  strain  at  the  inboard  (2.67  ft) 
station.  A  comparison  of  Figures  28  and  29  indicates  that  the  peak-to- 
peak  amplitude  of  the  longitudinal  strain  remained  essentially  the  same, 
but  the  frequency  content  was  appreciably  altered  at  both  stations. 
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Even  under  the  restrictive  assumption  noted  above,  the  experi¬ 
mental  strains  showed  appreciable  variations  in  magnitude  and  frequency 
content  over  the  time  interval  recorded.  In  some  cases,  the  peak-to- 
peak  amplitude  more  than  doubled  and  the  dominant  frequency  shifted 
from  IP  to  3P.  Because  of  the  nonperiodicity  of  the  strains  and  the 
sequential  method  of  collecting  the  strain  data,  noted  previously,  the 
resolved  longitudinal  and  shear  stress  values  are  highly  questionable, 
as  are  the  harmonics  of  the  longitudinal  strains. 
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LONGITUDINAL  STKCS5.  Is  tin  »  SH EAR  STRTSS,  Ik  /m  1  LONGITUDINAL  STRAIN.  in  /in  (atO1) 


RADIUS  =  2.67  ft 

MEAN  LONGITUOIRAL  STRESS  «  2017.66  Ik  /in 


RAOII*  =  3.67  (I 

MEAN  LONGlTUOlNAL  STRESS  ■  Jk27  36  Ik  /■«  1 


160 1  US  =  1.67  ft 

MEAN  LONGITUDINAL  STRAIN  0.3367  •  10'  *  in  /in 


V.  AZIMUTH  ANGLE,  dkg. 


f  .  AZIMUTH  ANGLE,  dtf. 


Figure  26.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  STRESSES  AND 
STRAINS  AT  TWO  RADIAL  STATIONS  FOR  =  39.8  ft  /sec  ; 

7  =  30';  n  =  1264  r.p.m. 
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Figure  27.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  STRESSES  AND 
STRAINS  AT  TWO  RADIAL  STATIONS  FOR  =  39.8  ft  / sec  ; 

T  =  40°;  fi  =  1264  r.p.m. 
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LONGITUDINAL  STRESS,  lb  /in  3  (,I0‘2)  SHEAR  STRESS,  lb  /in  3  LONGITUDINAL  STRAIN,  in  /in  («I03) 


Figure  28.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  STRESSES  AND 
STRAINS  AT  TWO  RADIAL  STATIONS  FOR  Vf  =  39.8  ft  /sec  ; 

T  =  50*;  /I  =  1264  r.p.m.;  WITHOUT  FLAPS  AND  SLATS 
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<¥  AZIMUTH  ANGIE.  deg. 


v>  .  AZIMUTH  ANGLE,  deg. 


Figure  29.  MEAN  AND  INSTANTANEOUS  VALUES  OF  EXPERIMENTAL  STRESSES  AND 
STRAINS  AT  TWO  RADIAL  STATIONS  FOR  Vf  =  39.8  ft  /sec  ; 
r  =  50  ;  n  =  1264  r. p.m. ;  WITH  FLAPS  AND  SLATS 


MEAN  LONGITUDINAL  STRAIN,  in  /in  (xIO1*) 


Figure  30.  EXPERIMENTAL  MEAN  LONGITUDINAL  STRAIN  AT  TWO  RADIAL 
STATIONS  VS.  ANGLE  OF  ATTACK  OF  THRUST  AXIS 
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